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ABSTRACT 


A  comprehensive  numerical  model  has  been  developed  to  investigate 
the  characteristics  of  an  impulsively-started,  asymptotically-steady, 
uniform  flow  about  a  circular  cylinder.  The  model  is  based  on  the 
rediscretization  of  the  shear  layers,  wake-boundary- layer  interaction, 
and  the  dissipation  of  vorticity  by  turbulence. 

The  forces  acting  on  the  cylinder,  rate  of  vorticity  flux, 
Strouhal  number,  oscillations  of  the  stagnation  and  separation  points, 
longitudinal  and  transverse  spacing  of  the  vortices,  and  the  base 
pressure  have  been  calculated  and  shown  to  be  in  conformity  with  those 
obtained  experimentally. 

The  model  has  been  used  to  predict  the  characteristics  of  hydro- 
elastic  oscillations  of  a  cylinder  in  the  range  of  synchronization. 
The  numerical  experiments  shed  considerable  light  on  the  interaction 
between  the  fluid  motion  in  the  wake  and  the  dynamics  of  the  body. 

An  extensive  sensitivity  analysis  has  been  carried  out  to  determine 
the  stability  of  all  the  parameters  and  hence  the  stability  of  the 
numerical  model  itself. 
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I.   INTRODUCTION 

A.   FLOW  SEPARATION 

Real  fluids  flowing  over  bluff  bodies  separate  under  the  influence  of 
adverse  pressure  gradient  as  momentum  is  consumed  by  wall  shear.  Separa- 
tion is  best  understood  in  terms  of  its  consequences.   In  fact,  it  was  the 
keen  observations  of  separation  in  d iff users  that  led  Prandtl  to  introduce 
the  concept  of  transition  or  boundary  layer  theory.   Even  though  the 
boundary  layer  theory  has  revolutionized  fluid  dynamics,  the  phenomenon 
which  gave  impetus  to  its  discovery  still  remains  an  enigma. 

Through  the  years  various  characteristics  of  separated  flow  about 
bluff  bodies  in  general  and  about  circular  cylinders  in  particular  have 
been  partially  understood  through  observations,  measurements,  and  analyses. 

The  separation  point  may  be  mobile  as  on  a  circular  cylinder  or   fixed 
along  a  line  as  on  a  wedge  or  blunt-based  body.   When  a  separation  point 
is  mobile,  the  pressure  decreases  up  to  the  point  of  maximum  velocity 
(slightly  upstream  of  separation)  and  then  increases  up  to  the  point  of 
separation.   When  the  separation  point  is  fixed,  the  pressure  continues 
to  decrease  up  to  the  point  of  separation.   Consequently,  there  are  signi- 
ficant differences  in  the  separation  of  flow  about  a  circular  cylinder 
and  a  sharp-edged  body.   In  either  case  the  velocity  distribution,  rate 
of  vorticity  flux,  and  the  amplitude  and  frequency  of  oscillation  of  the 
separation  points  are  not  independent  of  the  wake  in  the  downstream  side 
of  the  body.   It  is  this  interaction  that  makes  the  calculation  of  sepa- 
rated flows  exceedingly  complex.   In  other  words,  one  needs  the  evolution 
of  the  boundary  layer  on  the  forebody  (upstream  of  separation)  of  the 


cylinder  in  order  to  determine  the  instantaneous  position  of  the  separa- 
tion points.  But,  the  characteristics  of  the  boundary  layer  are,  in  turn, 
strongly  affected  by  the  fluid  motion  in  the  wake.   Evidently,  one  needs 
some  kind  of  an  interactive  mechanism  whereby  the  effect  of  the  wake  on 
the  boundary  layer,  and  hence  separation,  can  be  given  proper  considera- 
tion.  In  mathematical  terms  one  can  state  that  the  unseparated  potential 
flow,  or  the  d'Alembert  flow,  does  not  even  constitute  a  first  approxima- 
tion to  the  separated  real  flow.  Partly  the  realization  of  this  fact, 
partly  the  difficulties  associated  with  the  integration  of  the  Navier- 
Stokes  equations  across  a  singularity  (separation  point),  and  partly  lack 
of  sufficient  computational  capability  have  prevented  the  prediction  of 
the  kinematic  and  dynamic  characteristics  of  separated  flows  about  bluff 
bodies.   In  the  following,  the  basic  characteristics  of  flow  about  circu- 
lar cylinders  are  described  with  special  emphasis  on  those  aspects  which 
prepare  the  reader  for  the  understanding  of  the  numerical  model  developed 
during  the  course  of  this  work  and  the  assessment  of  the  similarities  and 
dissimilarities  between  the  predictions  and  observations  and  measurements. 

B.   VORTEX  SHEDDING 

The  separation  of  flow  gives  rise  to  alternate  vortex  shedding.   This 
phenomenon  is  intrinsic  to  the  flow  itself  and  has  been  known  at  least 
since  the  times  of  Leonardo  de  Vinci.   In  1878,  almost  exactly  100  years 
ago,  Strouha I  discovered  that  there  is  a  relationship  between  the  frequency 
of  vortex  shedding,  the  velocity  of  flow,  and  the  diameter  of  the  cylinder. 
This  relationship,  which  is  now  known  as  the  Strouhal  Number,  St,  is 
given  by  St  =  fD/U.   Strouhal 's  experiments  gave  an  average  value  of  St  - 
0.185.   During  the  course  of  numerous  experiments  since  then,  it  has  been 
experimentally  demonstrated  that:   (i)  the  Strouhal  number  is  dependent  on 


the  Reynolds  number  and,  as  first  suggested  by  Lord  Rayleigh,  must  be 
written  as  f  =  St(Re)U/D;  (ii)  the  vortex  shedding  frequency  is  only  an 
average  quantity  for  a  given  flow  and  that  there  are  always  secondary  and 
tertiary  frequencies  present;  and  finally,  (iii)  there  are  regions  of  flow 
such  as  the  region  of  drag  crisis  where  it  is  quite  difficult  to  identify 
a  particular  frequency.   It  is  evident  from  these  facts  that  in  some  re- 
gions of  flow  one  must  speak  of  the  spectral  content  of  the  wake  fluctua- 
tions in  assessing  the  value  of  the  Strouha I  number  and  the  forces  exerted 
on  the  cylinder.  The  foregoing  discussion  of  the  Strouha I  number  raises 
two  questions:  firstly,  is  it  possible  to  devise  a  universal  Strouha!  num- 
ber which  would  be  applicable  to  all  two  and  three-dimensional  bluff 
bodies,  and  secondly,  is  it  possible  to  theoretically  predict  the  Strouhal 
number  through,  for  example,  the  perturbation  of  the  d'Alembert  flow.  An 
attempt  along  these  lines  has  been  made  by  Sacksteder  (J  11  who  found  St  = 
0.2028  for  an  infinitely  large  Reynolds  number.   The  reasoning  behind 
this  prediction  and  its  validity  will  not  be  discussed  here  further. 

The  former  question  of  the  possibility  of  devising  a  universal  Strouhal 
number  has  been  of  interest  to  various  researchers.  Table  I,  shown  below, 
lists  those  definitions  together  with  their  originators. 

TABLE  I 


St  = 

Va  1  ue 

Body  Type 

Origi  nator 

f  h'/UK 

V 

0.  163 

two-dimensional  bodies 

Roshko  l2_1 

f  I   /u 

V  V 

0.181 

two-dimensional  bodies 

Bearman  [!>"] 

f  I   /K-U 
V  w 

0.190 

ax i symmetric  bodies 

Calvert  [4j 

f  I   /K-U 

V  w 

0.163 

two-dimensional  bodies 

Simmons  L5] 

In  Table  I,  h'  represents  the  wake  width  as  determined  from  Roshko' s 
notched-  hodograph  theory  E2~J;  I    ,    lateral  spacing  of  the  vortices  as 
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determined  by  Kronauner's  minimum  drag  criteria  as  reported  by  Bearman  [3]; 

I   ,  the  wake  width  at  a  distance  of  the  wake  establishment  region;  and  K  = 
w 


/P-CT  where  C  ,  is  the  base  pressure  coefficient.  All  these  definitions 
pb       pb 

require  either  the  measurement  or  the  solution  of  one  or  more  characteris- 
tics of  the  wake.  Thus,  their  capability  to  predict  a  vortex  shedding  fre- 
quency is  quite  limited. 

It  is  advantageous  to  note  at  an  early  stage  those  factors  which  in- 
fluence the  Strouhal  number.   Firstly,  small  amplitude  transverse  oscilla- 
tion of  cylinders  or  cables  distinctly  regularize  the  vortex  shedding  and 
hence  the  Strouhal  number.   For  example,  the  experiments  conducted  at  the 
National  Physical  Laboratory  [_6~]   have  shown  that  the  Strouhal  number  even 
in  the  drag  crisis  region  remains  perfectly  constant  and  equal  to  about 
0.21.   Secondly,  the  vortices  do  not  peel  off  from  a  cylinder  as  perfectly 
two-dimensional  line  vortices.  As  it  was  first  shown  by  Humphreys  L7_l, 
there  is  a  cell  pattern  separated  by  a  correlation  length.  The  longer  the 
cylinder  the  larger  the  number  of  cells.  The  consequences  of  the  forma- 
tion of  such  cells  are  that  there  are  phase  shifts  in  pressure  along  the 
cylinder  and  hence  a  reduction  in  the  integrated  transverse  force.   In 
other  words,  the  transverse  force  calculated  through  the  integration  of 
the  pressure  distribution  at  a  given  section  is  not  identical  to  the  aver- 
age transverse  force  determined  experimentally.   All  numerical  analyses 
which  model  the  complex  three-dimensional  flow  as  a  two-dimensional  flow 
must  necessarily  predict  a  transverse  force  or  lift  coefficient  somewhat 
larger  than  those  found  by  measurements.   Other  consequences  of  the  cell 
formation  concern  the  strength  of  the  vortices.   The  vortices  are  not 
necessarily  oriented  in  a  plane  normal  to  the  flow  and  the  measurements 
(flow  visualization  or  electro-optical  techniques)  yield  only  that 


component  of  vorticity  which  is  normal  to  the  plane  of  flow.  Additional- 
ly, experimental  results  necessarily  depend  on  the  length-to-diameter 
ratio,  end  conditions,  intensity  and  scale  of  turbulence  in  the  ambient 
flow,  the  rigidity  of  the  cylinder  mounting,  and  of  course  on  Reynolds 
number  and  the  surface  conditions  of  the  cylinder.  Most  recent  experi- 
ments (Shih  and  Hove  L"8_l)  have  shown  that  roughness  affects  the  base 
pressure  in  the  same  manner  as  the  small  transverse  oscillations  (beyond 
a  threshold  I  imit  of  about  0.10  diameters).   Uniform  roughness  yields  a 
perfectly  uniform  base  pressure  and  removes  the  cell  pattern. 

C.  RESISTANCE  IN  STEADY  FLOWS 
I  .   In-Line  Forces 

The  understanding  of  the  fluid  forces  exerted  on  bluff  bodies 
comes  primarily  from  measurements.   Of  all  bluff  bodies,  the  circular  cy- 
linder has  been  subjected  to  the  largest  number  of  investigations.  Those 
experiments  which  constitute  milestones  in  the  determination  of  the  magni- 
tude and  characteristics  of  resistance  are  the  works  by  Wiesel  sberger  L"9_l, 

Fage  and  Johansen  ClO],  and  Rhosko  Dl3.   A  plot  of  the  drag  coefficient, 

2 
Cn  defined  by  C  =  2F/pLDU  ,  is  shown  in  Fig.  I  as  a  function  of  the 

Reynolds  number  defined  by  Re  =  UD/v.   It  is  a  well  known  fact  that  for 

4  5 

Reynolds  numbers  larger  than  about  10  and  smaller  than  about  2x10  ,  the 

drag  coefficient  for  a  circular  cylinder  remains  essentially  constant  at 

5 
a  value  of  about  1.2.   For  Re  >  2x10  ,  the  shear  layers  emanating  from 

the  separation  points  become  turbulent.   Since  a  turbulent  flow  can  ex- 
change larger  momentum  between  the  inner  and  outer  flow  (i.e.,  between  the 
boundary  I ayer  and  externa  I  flow),  it  can  sustain  a  larger  adverse  pressure 
gradient  and  even  reattach  to  the  cylinder.  The  reattachment  of  the  tur- 
bulent shear  layers  gives  rise  to  the  formation  of  so  called  separation 
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bubbles  and  to  the  displacement  of  the  attached  flow  further  downstream 
from  the  first  separation  point.  The  reattached  flow  cannot  continue  to 
remain  attached  indefinitely  and  reseparates  at  an  angle  of  about  130  de- 
grees from  the  front  stagnation  point.  Thus,  the  transition  of  the  flow 
in  the  shear  layers  from  a  laminar  to  a  turbulent  state  results  in  a  con- 
siderably smaller  wake  and  hence  in  a  considerably  reduced  drag  coefficient. 
This  particular  flow  regime  is  commonly  called  the  drag  crisis  region  or 
the  critical  region.  As  the  Reynolds  number  increases,  parts  of  the  bound- 
ary layer  adjacent  to  the  first  separation  point  become  turbulent  and  the 
second  separation  point  begins  to  move  gradually  upstream.  This  region  is 
called  the  transcritica I  region.  At  higher  Reynolds  numbers,  the  entire 
boundary  layer  becomes  turbulent  and  the  drag  coefficient  reaches  a  nearly 
constant  value  of  about  0.6.  There  is  considerable  scatter  in  the  data 
in  the  supercritical  region  because  of  the  difficulties  associated  with 
flows  at  such  high  speeds. 

The  present  knowledge  of  the  drag  coefficient  for  a  smooth  cylin- 
der is  limited  to  Reynolds  numbers  smaller  than  about  10  .   There  is  great 
practical  interest  in  the  determination  of  force  coefficients,  i.e.,  drag 
coefficient,  lift  coefficient,  and  the  Strouhal  number  for  projects  such 
as  ocean  thermal  energy  conversion   (OTEC),  large  cooling  towers,  high  rise 
structures,  etc.  A  careful  consideration  of  the  various  possible  methods 
to  obtain  data  at  Reynolds  numbers  in  the  order  of  10  or  10   shows  that 
there  are  practically  insurmountable  difficulties.   The  existing  wind  and 
water  tunnels  are  not  capable  of  yielding  such  high  Reynolds  numbers  with- 
out the  adverse  effects  of  blockage  and  compressibility.   The  pop-up  and 
drop-down  experiments  with  large  cylinders  are  constrained  by  cavitation, 
weight,  and  stability  conditions.   It  might  be  possible  in  the  future  to 
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close  the  high  Reynolds  number  gap  by  building  large  cryogenic  wind  tun- 
nels operating  with  freon  or  liquid  helium. 

The  drag  coefficient  depicted  in  Fig.  I  may  be  influenced  by  a 
number  of  factors  such  as  roughness  and  turbulence.  Numerous  investiga- 
tors, notably  Fage  and  Warsap  D2_l,  Achenbach  Cl3l,  and  Guven  Cl4_l,  have 
clearly  demonstrated  that  the  larger  the  surface  roughness,  the  sooner  is 
the  occurrence  of  the  critical  region  and  the  higher  the  dip  in  the  drag 
coefficient  (see  Fig.  2).  Turbulence  plays  a  similar  role  in  the  evolu- 
tion of  the  drag  coefficients.   It  should  be  noted  from  Fig.  2  that  the 
drag  coefficient  for  a  rough  cylinder  at  supercritical  Reynolds  numbers  is 
about  twice  that  for  a  smooth  cylinder.   Furthermore,  the  experiments  by 
Guven  [14]]  and  Szechenyi  D  5_1  have  shown  that  at  sufficiently  high 
Reynolds  numbers  there  is  a  Reynolds  number  independent  region.   This  ex- 
perimental fact  has  been  termed  the  independence  principle  and  used  (Shih 
and  Hove  L~83)  to  extrapolate  the  drag  coefficient  data  to  Reynolds  numbers 

Q 

in  the  order  of  10  .   If  the  validity  of  the  independence  principle  can 
be  demonstrated  unambiguously,  it  may  obviate  the  need  for  difficult  ex- 
periments in  cryogenic  tunnels. 
2.  Transverse  Forces 

Lift  or  transverse-force  coefficient  data  are  far  less  numerous  and 
exhibit  a  great  deal  of  scatter.   Figure  3,  which  represents  most  of  the 

existing  lift  coefficient  data,  shows  that  the  lift  coefficient,  C.  defined 

2 
by  C  =  (amplitude  of  lift  force)/0.5pLDU  ,  may  vary  from  0.2  to  1.0  in  the 

subcritical  Reynolds  number  range.   In  the  critical  range  it  is  impossible 

to  define  a  lift  coefficient  without  resorting  to  its  root-mean-square  (rms) 

values.   It  is  for  this  reason  that  some  of  the  data  appearing  in  Fig.  3 

are  in  terms  of  the  rms  values  of  the  lift  coefficients.   For  a  more 
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detailed  discussion  of  the  spectral  content  of  lift  in  the  critical  region 
see  Bubl itz  [17]. 

The  reasons  associated  with  the  large  variations  of  the  lift  co- 
efficient are  directly  related  to  the  mobility  of  the  separation  points, 
formation  of  vortex  cells,  lack  of  spanwise  correlation,  end  effects,  tur- 
bulence present  in  the  ambient  flow,  etc.  Since  each  experiment  is  carried 
out  with  a  different  length-to-diameter  ratio,  blockage  ratio,  and  end 
conditions,  the  number  of  vortex  cells  over  the  cylinder  and  the  effect  of 
the  flow  between  the  end  of  the  cylinder  and  the  tunnel  wall  give  rise  to 
widely  different  lift  coefficients. 

The  effect  of  roughness  on  the  lift  coefficient  is  to  reduce  its 
scatter  and  to  increase  its  magnitude.  As  noted  earlier,  roughness  dis- 
tinctly regularizes  the  vortex  formation  and  the  base  pressure  thereby 
eliminating  the  effect  of  phase  shifts  which  cause  the  said  scatter  in 
lift.  Data  supporting  the  above  experimental  facts  will  not  be  reported 
here  for  sake  of  brevity.  For  additional  details  the  reader  is  referred 
to  the  work  by  Szechenyi  Cl53. 

D.   RESISTANCE  IN  TIME-DEPENDENT  FLOWS 

There  are  an  infinite  number  of  time  dependent  flows  and  it  is  impos- 
sible to  discuss  their  characteristics  in  any  unifying  manner.  One  might 
obtain  a  clearer  insight  into  such  flows  by  considering  only  those  cases 
which  are  relatively  more  manageable  for  numerical  or  laboratory  experiments, 
Among  such  flows  which  have  been  subjected  to  extensive  numerical  analysis 
and  experiments  are  the  impulsively  started  flows,  uniformly  accelerating 
■flows,  and  harmonically  oscillating  flows.   Obviously,  no  distinction  is 
made  here  between  the  said  fluid  flows  about  a  body  at  rest  and  a  body 
undergoing  similar  time  dependent  motions  in  a  fluid  otherwise  at  rest. 
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I  .   Impulsively  Started  Flows 

The  impulsively  started  flow  has  been  by  far  the  one  most  investi- 
gated numerically  through  the  use  of  the  Navier-Stokes  equations.   These 
investigations  (Collins  and  Dennis  Cl8H,  and  Yang  Cl9j)  have  been  confined 
to  the  very  early  stages  of  the  motion  during  which  the  separation  points 
move  from  the  rear  stagnation  point  to  about  109  degrees.  An  exception 
to  the  foregoing  is  the  work  reported  by  Thoman  and  Szewczyk  C20H.  They 
have  solved  the  Navier-Stokes  equations  through  the  use  of  a  suitable 
finite  difference  technique  for  various  Reynolds  numbers  and  found  that 
the  drag  coefficient  overshoots  for  a  relative  flow  displacement  (normalized 
time)  of  Ut/c  =  5  and  then  gradually  reaches  steady  state  values. 

Only  two  separate  investigations  of  the  impulsively  started  flow 
have  been  reported.  The  first  is  that  carried  out  by  Schwabe  C2lJ  upon 
Prandtl's  suggestion  and  the  other  by  Sarpkaya  E22,  23].   Schwabe's  experi- 
ments were  conducted  with  a  9  cm  diameter  cylinder  at  a  Reynolds  number 
of  about  600.   His  drag  coefficient,  calculated  indirectly  through  the  use 
of  motion  pictures,  rapidly  reaches  a  value  of  about  2.0  and  remains  con- 
stant thereafter.   Schwabe's  findings  are  not  in  conformity  with  either 
the  numerical  calculation  or  with  the  experimental  results  obtained  by 
Sarpkaya.   Figure  4  shows  three  sets  of  representative  data  at  Reynolds 
numbers  from  16000  to  35000  for  the  drag  coefficient.  Clearly,  the  very 
initial  instants  of  the  motion  cannot  be  analyzed  because  of  the  problems 
associated  with  the  generation  of  an  impulsively  started  flow.   At  Ut/c 
slightly  greater  than  zero,  the  force  acting  on  the  cylinder  is  almost 
entirely  inertial.   As  the  separation  sets  in  at  the  rear  stagnation  point 
and  rapidly  moves  upstream,  the  drag  also  increases  rapidly  and  Cn  reaches 
a  value  of  3bout  1,5  at  Ut/c  between  4.0  and  5.0.   Then  C  drops  rapidly 
to  the  steady  state  value  for  smooth  cylinders  in  steady  uniform  flows  at 
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Ut/c  larger  than  about  7.0.  Detailed  flow  visualization  experiments  by 
Sarpkaya  [22]  have  shown  that  the  reason  for  the  overshoot  in  C„  is  the 
symmetric  development  of  two  vortices  in  the  cylinder  wake  and  the  rapid 
accumulation  of  vorticity  in  each  vortex.   It  must  be  noted  that  for  such 
a  symmetric  vortex  configuration,  the  only  means  by  which  the  vorticity 
may  be  cancelled  are  the  exchange  of  vorticity  between  the  two  vortices 
and  the  cancellation  of  vorticity  with  oppositely-signed  vorticity,  so 
called  counter-vorticity,  generated  in  the  rear  boundary  layers.   In  other 
words,  the  rate  of  accumulation  of  vorticity  in  the  vortices  is  consider- 
ably larger  than  that  at  any  other  time.  Other  vorticity  cancellation 
mechanisms  which  come  into  play  during  the  separation  of  the  vortex  from 
its  feeding  sheet,  as  proposed  by  Gerrard  [24],  do  not  play  any  role  dur- 
ing the  growth  of  the  symmetric  vortices.  Thus,  the  rapid  accumulation 
of  vorticity  gives  rise  to  the  drag  overshoot.   As  the  asymmetry  sets  in, 
vortices  shed  a  I ternati ngly  and  the  flow  gradually  becomes  asymptotically 
steady.  The  possible  mechanisms  which  cause  the  asymmetry  in  an  otherwise 
symmetrical  flow  about  a  symmetrical  body  are  discussed  in  greater  detail 
later.   Suffice  it  to  note  that  one  of  the  fundamental  facts  of  fluid 
flows  is  that  symmetric  causes  do  not  give  rise  to  symmetric  effects,  as 
most  eloquently  philosophized  by  Birkhoff  [253. 

Following  the  onset  of  asymmetry,  the  cylinder  experiences  a 
transverse  force,  as  shown  in  Fig.  5,  in  terms  of  the  lift  coefficient  and 
the  normalized  time  Ut/c.   The  lift  force  fluctuates  with  the  Strouha I  fre- 
quency.  The  shedding  of  each  vortex  gives  rise  to  a  small  fluctuation  in 
CL.  The  frequency  of  this  fluctuation  is  twice  the  Strouhal  frequency 
since  two  vortices  (one  from  top  and  from  bottom)  are  shed  during  a  given 
cycle  of  I  ift  osci  I lations. 
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2.  Other  Time-Dependent  Flows 

As  noted  earlier,  the  uniformly  accelerating  flows  and  harmonical- 
ly oscillating  flows  have  been  subjected  to  numerous  analytical  and  experi- 
mental investigations.  Among  those  notable  are  the  investigations  carried 
out  by  Iverson  and  Balent  C263,  Keim  [27],  Hamilton  C283,  Odar  and  Hamilton 
C29H,  Sarpkaya  and  Garrison.  C30H,  and  Sarpkaya  C3l].  The  basic  question 
has  been  and  continues  to  be  the  formulation  of  a  resistance  equation  which 
takes  into  consideration  the  effects  of  acceleration,  velocity,  and  the 
history  of  the  motion.  Ail  investigations  have  shown  that  resistance  in 
time-dependent  flows  in  general  and  in  harmonically  oscillating  flows  in 
particular  cannot  be  considered  as  a  juxtaposition  of  resistance  in  steady 
flows  at  corresponding  instantaneous  Reynolds  numbers.  Stokes,  in  his 
celebrated  paper  on  the  sinusoidal  oscillation  of  a  sphere  in  a  fluid 
otherwise  at  rest,  has  shown  that  the  fluid  resistance  is  comprised  of  a 
linear  velocity-dependent  force  and  an  acceleration-dependent  inertial 

force.  The  coefficients  for  both  forces  depend  on  the  Reynolds  number  and 

2 
the  frequency  parameter  D  /vT  where  T  is  the  period  of  oscillation.  Ob- 
viously, Stokes'  analysis  is  valid  for  only  very  small,  as  yet  unspecified, 
Reynolds  numbers.   Inspired  partly  by  Stokes'  contribution  and  partly  by 
practical  considerations,  Morison  L32j  proposed  a  resistance  equation  which 
may  be  written  as 

F=   2"pAp   CDUIUI    +PVCMdT  (l) 

in  which  U  represents  the  instantaneous  velocity;  A  ,  the  projected  area; 
V,  the  displaced  volume;  and  Cn  and  CM,  the  drag  and  inertia  coefficients 
respectively.   For  harmonically  oscillating  flow  about  a  stationary  cylin- 
der where  the  fluid  velocity  is  U  =  -11.  cos  8  (here  8  =  2lTt/T)  Eq.  (I) 


30 


reduces  to 


1  =  -  Cjcos   el  cos   8  +  ;,  T/n'C,.  sin   9  (2) 


ipDLUM2      D'      ' 

The  essence  of  this  equation  is  that  the  resistance  in  harmonic  flow  is 
assumed  to  be  equal  to  a  linear  combination  of  a  velocity-squared  dependent 
drag  force  and  an  acceleration-dependent  inertia  force.   In  the  midst  of 
insufficiently  clear  insight  into  the  kinematics  of  the  complex  flow,  it 
has  not  yet  been  possible  to  demonstrate  the  reasons  for  the  differences 
between  the  measured  force  and  that  calculated  from  Morison's  equation  for 
ILT/D  values  in  the  vicinity  of  15.  Additional  experiments  and  numerical 
analysis,  possibly  through  the  use  of  the  discrete  vortex  model,  might  shed 
further  light  on  the  subject  and  into  the  overall  understanding  of  resis- 
tance in  time-dependent  flows.   It  is  evident  from  the  discussion  of  both 
the  steady  and  unsteady  flows  that  analytical  or  numerical  methods  are 
urgently  needed  with  which  relatively  inexpensive  investigations  can  be 
conducted  to  clarify  the  complex  wake-boundary- I ayer  interaction  mechanism. 
Of  all  the  possible  flows  which  might  be  considered  as  a  candidate  for  such 
an  analysis,  the  impulsively  started  flow  is  certainly  the  one  most  manage- 
able and  most  fundamental  to  the  understanding  of  ail  other  flows. 

E.  METHODS  OF  NUMERICAL  ANALYSIS 

With  the  advent  of  high  speed  computers,  numerical  methods  began  to 
play  a  greater  role  in  fluid  dynamics.   Currently  four  major  methods  of 
analysis  are  in  use:   finite  difference,  finite  element,  marker-cell  tech- 
nique, and  the  discrete  vortex  method. 

Finite  difference  methods  have  been  used  for  the  solution  of  flow 
about  a  circular  cyl  inder  for  relatively  small  Reynolds  numbers  (Re  -    100) 
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either  through  the  use  of  the  Navier-Stokes  equations  as  they  are  commonly 
written  or  through  the  more  frequently  used  vorticity  equations.   Ingham 
L~33H  calculated  the  steady  flow  past  a  circular  cylinder  at  Re  =  100. 
Thoman  and  Szewczyk  E20H,  as  noted  earlier,  carried  out  a  similar  numeri- 
cal analysis  of  flow  over  stationary  and  rotating  cylinders  at  Reynolds 
numbers  as  large  as  10  .  There  are  a  number  of  such  examples  in  the 
literature  but  they  are  not  cited  here  because  of  their  lack  of  relevance 
to  the  subject  under  consideration. 

The  finite  element  method,  which  came  into  prominence  first  in  the 
area  of  structural  analysis  and  then  permeated  through  heat  transfer  and 
fluid  mechanics,  has  been  primarily  used  for  unseparated  flows  about 
bluff  bodies.  As  far  as  fluid  mechanics  is  concerned,  the  outstanding  ex- 
amples of  such  applications  may  be  found  in  two  recent  volumes  of  work 
dealing  with  finite  elements  in  fluid  mechanics  L~34_l.   Bretanow  and  Ecer 
L~35j  made  an  attempt  to  calculate  the  separated  flow  about  a  circular 
cylinder  at  small  Reynolds  numbers  but  the  results  have  been  inconclusive. 
Applications  of  the  finite  element  method  to  flows  about  foils  where  the 
flow  remains  attached  have  yielded  results  in  conformity  with  experiments 
(Bretanow  and  Ecer  C35H). 

The  marker-cell  technique  developed  by  the  researchers  at  the  Los 
Alamos  Scientific  Laboratory  during  the  past  20  years  uses  the  finite 
difference  form  of  the  Navier-Stokes  equations  and  then  adjusts  the  pres- 
sure in  each  cell  to  insure  that  the  equation  of  continuity  is  satisfied. 
Numerous  ingenious  techniques  have  been  developed  by  Fromm  and  Harlow 
[36]  and  Nichols  and  Hirt  L"37_]  for  the  efficient  adjustment  of  pressure 
in  each  cell.   Such  techniques  will  not  be  dealt  with  here.   The  researchers 
of  Los  Alamos  Scientific  Laboratory  have  applied  the  marker-cell  technique 
to  the  flow  past  cylinders,  spheres,  rectangular  blocks  and  to  various  other 
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free-surface  flows.  These  applications  dealt  with  very  low  Reynolds  num- 
bers (in  the  order  of  100).  Even  though  the  technique  has  considerably 
increased  the  understanding  of  the  kinematics  of  the  wake,  the  precise 
nature  of  the  time-dependent  force  acting  on  any  one  of  the  bodies  cited 
above  remains  obscure. 

The  discrete  vortex  model,  which  is  included  as  one  of  the  four  major 
models,  is  the  subject  of  extensive  discussion  in  the  next  chapter.  There- 
fore, the  remainder  of  this  section  is  confined  to  a  discussion  of  the  suc- 
cesses and  shortcomings  of  the  first  three  methods. 

All  numerical  models  are  by  their  very  nature  approximate.  They  satis- 
fy certain  critical  conditions  only  at  discrete  time  steps.  Secondly,  they 
require  mesh  or  element  sizes  which  vary  from  one  region  of  flow  to  another 
to  insure  that  the  characteristic  gradients  of  the  flow  (such  as  velocity 
or  pressure)  do  not  significantly  or  abruptly  change  over  any  one  element. 
For  example,  the  finite  difference  method  when  applied  to  the  viscous  flow 
over  a  circular  cylinder  requires  that  the  size  of  the  elements  on  and  near 
the  cylinder  be  considerably  smaller  than  the  boundary  layer  thickness. 
Otherwise,  the  evolution  of  the  boundary  layer,  the  positions  of  the  sepa- 
ration points,  and  the  evolution  of  the  wake  cannot  be  accurately  calcu- 
lated.  Since  the  boundary  layer  thickness  for  laminar  flow  is  proportional 

-1/2 
to  Re    ,  one  realizes  that  only  for  relatively  small  Reynolds  numbers  can 

a  sufficiently  large  boundary  layer  thickness  and  a  correspondingly  large 
mesh  size  be  obtained.   This  fundamental  difficulty  has  prevented  the 
finite  difference,  finite  element,  and  the  marker-cell  technique  from  making 
significant  progress  into  the  calculation  of  steady  or  time  dependent  flows 
about  bluff  bodies  at  Reynolds  numbers  of  practical  interest.   None  of  the 
three  techniques  cited  above  has  been  applied  to  the  solution  of  flow  prob- 
lems associated  with  oscillating  bodies.   The  reasons  for  this  are  rather 
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self-evident  since  the  solution  requires  not  only  a  moving  grid  but  also 
a  grid  whose  mesh  size  must  be  varied  with  time  to  minimize  the  effect  of 
sharp  velocity  and  pressure  gradients.   Finally,  it  should  be  noted  that 
all  three  numerical  techniques  present  "stiffer"  solutions  and  introduce 
an  artificial  numerical  viscosity  whose  influence  cannot  be  readily 
assessed.  Not  to  be  entirely  critical  of  the  said  methods,  for  there  are 
corresponding  shortcomings  of  the  discrete  vortex  model,  one  must  note 
that  in  the  wide  spectrum  of  Reynolds  numbers  there  will  always  be  some 
interest  in  low  Reynolds  number  flows.   Thus,  such  numerical  techniques 
will  continue  to  play  some  role  in  separated   flows  where  the  discrete 
vortex  model  is  not  advantageously  suited.   It  therefore  appears  that  no 
one  model  wi I  I  stand  out  as  superior  and  universal ly  appl  icable  and  that 
all  methods  will  continue  to  complement  each  other. 
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I  I .   FUNDAMENTALS  OF  THE  DISCRETE  VORTEX  MODEL 

A.   INTRODUCTION 

The  discrete  vortex  model  (DVM)  is  a  potential  flow  representation  of 
the  shear  layers  emanating  from  the  separation  points.  Any  separated  vis- 
cous flow  is  completely  specified  if  the  spatial  and  temporal  distribution 
of  vorticity  in  the  entire  flow  field  is  known.   For  reasons  described  in 
the  previous  chapter,  it  has  not  yet  been  possible  to  achieve  this  objec- 
tive through  any  exact  or  approximate  methods.   It  has  been  observed  and 
experimentally  verified  by  Fage  and  Johansen  DO]  that  the  shear  layers 
at  sufficiently  high  Reynolds  numbers  are  quite  thin  and  the  vorticity  is 
confined  into  regions  of  tightly  spiralled  vortex  sheets.   This  fact  alone 
invites  one  to  divide  the  shear  layer  into  a  number  of  small  segments  and 
concentrate  the  vorticity  of  each  segment  into  a  line  vortex  situated  at 
the  geometric  center  of  the  segment.   This  idealization  permits  one  to 
use  the  powerful  complex  function  theory  to  determine  the  kinematics  and 
dynamics  of  the  flow.   Evidently,  the  validity  of  the  results  can  only 
be  judged  in  light  of  experimental  facts. 

In  the  following,  only  the  gross  features  of  the  method  are  described 
and  the  details  are  presented  following  an  historical  resume. 

At  a  given  time  the  separated  shear  layers  are  represented  by  vortices 
situated  at  the  center  of  each  segment.  The  potential  flow  theory  permits 
one  to  calculate  the  velocity  of  each  vortex  induced  at  its  center  by  all 
other  vortices  and  the  ambient  flow.   Then  the  vortices  are  convected  for 
a  short  time  interval  At  through  the  use  of  a  suitable  convection  scheme 
to  new  positions.   It  is  then  necessary  to  introduce  new  vorticity  into 
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the  flow  at  or  near  the  separation  points  to  account  for  the  vorticity 
generated  in  the  boundary  layer.  The  introduction  of  new  vortices  (nascent 
vortices)  represents  an  addition  to  each  set  of  vortex  arrays.  Then  the 
process  is  repeated. 

In  spite  of  its  apparent  simplicity,  there  are  a  number  of  problems 
which  must  be  resolved  in  carrying  out  the  said  calculations:  (i)  the  vor- 
ticity introduced  at  each  time  step  must  be  related  to  the  characteristics 
of  the  boundary  layers  upstream  of  the  separation  points;  (ii)  one  must 
insure  that  the  point  vortices  continue  to  represent  the  vorticity  at  each 
segment  in  accordance  with  the  fundamentals  of  hydrodynamics;  (iii)  allow- 
ance must  be  made  for  the  interaction  of  oppositely-signed  vorticity  and 
the  decay  of  vortices;  (iv)  methods  must  be  devised  to  limit  the  number  of 
the  point  vortices  to  minimize  the  computation  time  without  affecting  the 
kinematics  and  dynamics  of  the  flow;  and  (v)  procedures  must  be  established 
for  the  flow  to  retain  its  inherent  flexibility  without  introducing 
artificial  stiffening  effects  (e.g.,  to  allow  for  the  mobility  of  stagna- 
tion and  separation  points). 

The  foregoing  conveys  the  idea  that  the  discrete  vortex  model  is  and 
will  continue  to  rely  on  a  strong  interaction  between  potential  and  vis- 
cous flow  theories.   Neither  is  capable  of  providing  a  solution  alone  but 
together  they  can  make  inroads  towards  describing  the  basic  features  of 
the  separated  flow  provided  that  sufficient  ingenuity  is  incorporated  into 
the  said  interaction.   Lastly,  it  must  be  noted  that  all  the  numerical 
methods  described  above  depend  on  the  availability  of  a  high-speed  digital 
computer. 

B.   HISTORICAL  DEVELOPMENT 

The  DVM  has  been  in  existence  for  the  past  46  years  since  its  incep- 
tion by  Rosenhead  [38]  in  1932.   Rosenhead  investigated  the  so-called 
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"Helmholtz  instability"  of  a  surface  of  discontinuity,  separating  two 
fluids  of  equal  speed  and  density  but  moving  in  opposite  directions.  His 
primary  objective  was  to  give  the  shear  layer  between  the  two  fluids  an 
initial  disturbance  by  situating  the  shear  layer  along  the  curve  y  =  0. I 
sin  (2irx/A)  and  then  observe  the  evolution  of  the  sheet  with  time.  To  ac- 
complish this,  Rosenhead  represented  the  shear  layer  by  point  vortices 
placed  along  the  sheet  at  intervals  such  that  their  projections  along  the 
x-axis  were  equ i-distant.   Each  vortex  was  assigned  a  strength  of  UA/6n 
where  n  is  the  number  of  vortices  per  wavelength, X,  and  U  the-  speed  of 
the  two  fluids.   Keeping  the  vortex  strengths  constant,  Rosenhead  followed 
the  deformation  of  the  sheet  by  convecting  the  vortices  with  a  simple 
Eulerian  integration  of  ri+A+  =  H  +  q^At  where  r^  is  the  position  vector 
of  the  j-th  vortex;  q^ ,  the  computed  velocity  vector;  and  At,  the  time 
interval.   Using  n  =  2,  4,  8,  and  12,  Rosenhead  observed  that  in  each  case 
the  sheet  rolled  up  smoothly  into  concentrated  clusters  of  vortices  posi- 
tioned at  intervals  of  X   along  the  x-axis.   Rosenhead's  results,  which  were 
computed  using  a  desktop  calculator,  are  shown  in  Fig.  6.  The  tendency  of 
a  distribution  of  point  vortices  to  roll  up  in  a  manner  consistent  with 
observations  is  clearly  evident. 

Rosenhead's  results  remained  undisputed  until  1959  when  the  computa- 
tional power  of  the  digital  computer  was  applied  to  the  analysis.  Birkhoff 
and  Fisher  [393  repeated  Rosenhead's  analysis  using  a  Runge-Kutta  integra- 
tion scheme,  22  vortices  per  wavelength,  and  smaller  time  intervals.   In- 
stead of  improved  results  they  observed  that  the  vortices  took  irregular 
paths  and  no  longer  rolled  up  smoothly.   Birkhoff  and  Fisher  concluded  that, 
given  sufficient  time,  an  originally  organized  array  of  equ i-strength  point 
vortices  would  eventually  become  random  and  that  the  validity  of  represent- 
ing a  shear  layer  by  point  vortices  is  questionable. 
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Fig.  6.   Example  of  Rosenhead 
calculations  f_  3  3  H  . 
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Shortly  thereafter,  Hama  and  Burke  C40H  undertook  a  similar  computer 
study  and  verified  the  observations  of  Birkhoff  and  Fisher.  However,  Hama 
and  Burke  noted  a  conceptual  error  in  Rosenhead's  original  method  of  plac- 
ing the  discrete  vortices  along  the  shear  layer.   Specifically,  if  the 
initial  vortex  sheet  was  to  have  a  constant  vorticity  per  unit  distance 
along  its  length  and  be  represented  by  equ i-strength  vortices,  then  the 
vortices  should  be  placed  at  equal  distances  measured  along  the  sheet  vice 
at  equal  x  increments.  Using  this  arrangement  Hama  and  Burke  repeated 
their  calculations  and  found  that  the  onset  of  irregular  vortex  motion  was 
significantly  delayed.   Eventually,  however,  a  random  pattern  did  result. 
Thus,  at  this  point  it  was  hypothesized  that  the  onset  of  instability 
could  be  delayed  although  the  questions  of  why  it  occurred  and  how  it 
could  be  controlled  were  still  a  mystery. 

Without  fully  realizing  the  consequences  of  sheet  instability  (the 
result  of  which  is  discussed  later),  the  DVM  continued  in  its  evolution 
relying  on  the  assumption  that  meaningful  results  could  be  obtained  up  to 
the  onset  of  irregular  vortex  motions.   In  1962  Abernathy  and  Kronauer  C4ll 
applied  the  method  to  investigate  the  motion  of  two  shear  layers  of  oppo- 
sitely-signed vorticity  separated  by  a  distance  h  and  each  given  an 
initial  perturbation  as  in  Rosenhead's  work.   They  observed  that  the  ver- 
tices accumulated  into  clusters  resembling  the  pattern  of  a  vortex  street. 
An  example  of  their  results  is  shown  in  Fig.  7. 

Encouraged  by  the  demonstrated  ability  of  the  DVM  to  represent  a  vor- 
tex street,  the  next  logical  extension  of  the  method  was  to  analyze  the 
flow  past  a  bluff  body.   However,  the  inclusion  of  a  solid  body  in  the 
flow  field  necessitated  new  considerations.   Specifically,  the  previous 
analyses  had  begun  with  a  discretized  sheet  having  a  preassigned  initial 
position  and  vorticity  distribution  along  its  length  (usually  taken  as 
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Fig.  7.   Example  of  Abernathy 

and  Kronauer's  calculations 
(h/X=0.34)  [41]. 
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constant).  The  application  of  the  DVM  to  flow  past  a  bluff  body  required 
that  the  discretized  shear  layers  be  generated  continuously  by  adding 
additional  point  vortices  into  the  flow  at  regular  time  intervals  and  at 
locations  representative  of  where  the  shear  layers  actually  emanate  from 
the  body.  The  requirement  for  some  method  to  account  for  vorticity  genera- 
tion and  introduction  in  an  otherwise  potential  flow  model  posed  a  signi- 
ficant problem.  Compounding  the  difficulty  was  the  time  dependent  nature 
of  both  the  generation  rate  and  its  point  of  introduction. 

The  first  DVM  analysis  to  include  a  bluff  body  was  performed  by 
Gerrard  G42]  for  a  circular  cylinder.  Although  a  circular  cylinder  was  of 
great  practical  interest  it  was  also  the  most  difficult  one  to  analyze 
because  of  the  unsteady  vorticity  generation  and  mobile  separation  points. 
Gerrard  was  forced  to  adopt  a  rather  arbitrary  method  of  introducing  new 
vortices  which  severely  influenced  the  results  of  his  analysis,  although 
it  still  represented  a  noteworthy  application  of  the  DVM.   Noting  the  dif- 
ficulties associated  with  the  mobile  separation  points  and  the  vorticity 
introduction,  Sarpkaya  [43]  introduced  a  mobile  separation  region  repre- 
sented by  an  infinite  number  of  vortex  sheets  of  vanishingly  small  vorticity 
connected  to  the  nascent  vortex.   This  analysis  was  restricted  to  impul- 
sively started  symmetric  flow  about  a  circular  cylinder.   In  the  ensuing 
years  many  investigators  bypassed  the  difficulty  introduced  by  the  mobility 
of  the  separation  points  and  dealt  with  bodies  with  sharp  edges  (e.g., 
inclined  plates,  blunt-based  bodies,  wedges,  etc.)  where  the  separation 
points  are  fixed.   Despite  the  growing  number  of  applications  of  the  DVM, 
a  number  of  disturbing  difficulties  remained  unresolved  and  a  full  inter- 
action between  the  boundary  layer,  separation  points,  and  the  wake  had 
not  been  established.   These  difficulties  required  the  use  of  a  number  of 
nondisposable  parameters.   Consequently,  the  proper  balance  between 
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heuristfc  reasoning  and  computational  formulation  had  to  evolve  through 
the  recognition  of  the  reasons  underlying  the  basic  difficulties. 

Briefly,  the  most  important  of  these  difficulties  stemmed  from  the 
proximity  of  the  vortices  to  each  other  and  to  the  solid  boundary,  and 
from  the  necessity  of  representing  the  continuous  process  of  vorticity 
generation  at  the  singular  separation  points  by  a  discrete  process.  The 
occasional  proximity  of  two  point  vortices  resulted  in  large  mutually- 
induced  velocities  and  hence  the  shooting  of  the  two  vortices  rapidly 
away  from  the  flow  field.  This  was  particularly  true  when  the  vortices 
ceased  to  be  located  at  the  geometric  center  of  the  sheet  segment  which 
they  were  supposed  to  represent.   Furthermore,  the  spiralling  of  the  sheets 
brought  vortices  on  two  adjacent  spirals  close  to  each  other  and  gave  rise 
to  a  traveling  and  growing  instability.  Thus,  it  was  not  possible  to 
correctly  represent  the  core  of  tightly  spiralled  sheets.  The  proximity 
of  the  vortices  to  the  boundary  and  hence  to  their  images  inside  the 
body  resulted  not  only  in  similar  instabilities  but  also  in  unrealistic 
local  pressures.  This,  in  turn,  resulted  in  unnatural  fluctuations  in 
drag  and  lift  forces. 

Faced  with  the  problems  just  cited  and  the  desire  to  enhance  the 
applicability  of  the  method,  various  investigators  devised  several  numeri- 
cal cures.   Chorin  and  Bernard  [441  introduced  a  cut-off  distance  from 
the  center  of  each  vortex  beyond  which  the  vortices  were  to  behave  like 
potential  vortices.   For  radial  distances  smaller  than  the  cut-off  distance 
the  vortices  were  assumed  to  exhibit  a  solid  body  rotation.   Their  calcu- 
lations did  not,  however,  result  in  an  improvement  of  the  analysis  of 
separated  flow  about  a  circular  cylinder.   Moore  [45],  in  dealing  with 
two  vortex  sheets  emanating  from  a  plane  lifting  surface,  introduced  the 
method  of  amalgamation.  When  a  vortex  sheet  was  spiralled  tightly  enough, 
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the  last  vortex  on  the  sheet  was  connected  to  a  core  vortex.  At  each  time 
step  the  last  vortex  was  amalgamated  into  the  core  vortex  increasing  its 
strength.  This  procedure  partly  avoided  the  difficulties  stemming  from 
the  proximity  problems  of  the  adjacent  coils  of  the  spiralling  sheet. 
Both  investigations  dealt  with  the  devising  of  methods  to  overcome  the 
instability  problems  rather  than  with  the  understanding  of  the  real  cause 
of  the  di  f f i  cu I ty. 

The  basic  question  of  why  an  initially  equidistant  set  of  vortices 
should  after  several  convections  come  in  close  proximity  of  each  other  was 
resolved  by  Fink  and  Soh  C46H  who  demonstrated  that  the  primary  source  of 
the  difficulty  was  in  failing  to  rediscretize  the  sheet  at  each  time  in- 
terval of  the  calculation.   In  other  words,  it  is  necessary  to  insure  that 
the  point  vortices  are  always  located  at  the  geometric  center  of  the  seg- 
ment of  the  sheet  which  they  are  supposed  to  represent.   The  complex  con- 
jugate velocity  q(z.)  at  a  point  z.  due  to  a  single  vortex  of  strength  I\ 


k 


located  at  z,  is  given  by 


q(z.)  =  J (3) 

M   J      2 1 TT  Z  .  -  Z. 

J        J    k 
For  a  series  of  N  discrete  vortices  on  a  vortex  sheet,  the  same  velocity 

at  a  point  z.  becomes  l~46H 
J 


N     F,       ,    F,  e 


J     2l  TT  L*l       Z.  -  Z,     2lTT   S.,./0-S. 


in 


k/j   j    k       i  j  +  l/2  J--I/2' 


zrzj+i/2 

Yzj-i/2 


(4) 


Mathematical  details  of  the  relevant  potential  flow  theory  will  be 
presented  in  Chp.  III.   Here  no  attempt  is  made  to  introduce  the  precise 
definitions  of  the  variables  appearing. 
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Equation  (4)  resembles  Eq.  (3)  for  the  induced  velocity  field  due  to  point 
vortices  except  for  the  logarithmic  term.  The  point  z.  lies  within  the 
segment  (s.  , ,~,  s.  ,,~)  without  necessarily  bisecting  it  at  all  times,  and 
s  measures  the  distance  along  the  sheet.  The  consequences  of  Eq.  (4)  are 
that:  (i)  if  the  equivalent  point  vortex  is  not  placed  at  the  midpoint  of 
its  segment  through  red iscretization  of  the  sheet  at  each  time  interval, 
then  the  logarithmic  term  does  not  vanish  and  the  computational  error  in- 
creases depending  on  the  number  of  vortices,  the  time  interval  used,  and 
the  total  time  of  computation;  (ii)  the  vortices,  which  initially  bisect 
the  segment  which  they  are  to  represent,  do  not  continue  to  do  so  at  the 
succeeding  time  intervals;  (iii)  the  use  of  finite  vortex  cores,  amalgama- 
tion of  vortices  at  the  center  of  the  spiral,  or  other  techniques  only 
delay  or  minimize  the  accumulation  of  the  errors  resulting  from  the  loga- 
rithmic term  in  an  amount  related  to  the  distance  between  z.  and  the  center 

J 

of  the  segment;  and  that  (iv)  the  growth  of  the  computational  error  may  be 

significantly  reduced  by  placing  each  discrete  vortex  at  the  midpoint  of 

its  segment,  i.e.,  by  p laci  ng  the  vortex  at  z.  =  0.5(z .  , ,  ,„+z  .  .  /0)  at  each 

J        J+l/2  j- 1/2 

time  interval.   Only  through  such  a  procedure  can  one  make  the  logarithmic 
term  vanish. 

The  calculations  are  then  carried  out  at  each  time  step  by  represent- 
ing the  vorticity  density  by  an  entirely  new  set  of  equidistant  vortices 
whose  strengths  are  adjusted  to  give  a  good  representation  of  that  density. 

The  foregoing  procedure  does  not  resolve  all  of  the  computational 
errors  particularly  in  regions  where  the  radius  of  curvature  of  the  sheet 
is  small.   Fink  and  Soh  applied  this  technique  to  the  early  stages  of 
relatively  simple  flows  with  fixed  separation  points.   Their  examples  did 
not  deal  with  the  alternate  shedding  of  vortices. 
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Even  though  the  method  of  rediscretization  essentially  resolved  the 
problems  stemming  from  the  proximity  effects,  two  fundamental  links  between 
the  potential  flow  model  and  the  behavior  of  real  fluids  remained  unre- 
solved. The  first  of  these  is  the  relationship  between  the  vorticity  to 
be  introduced  at  each  time  step  and  the  boundary  layer  which  generates  it. 
The  second  is  the  decay  of  vortices  by  various  mechanisms.  The  first 
question  was  initially  considered  by  Deffenbaugh  and  Marshall  f_47]  who  in- 
vestigated the  impulsively  started  flow  about  a  circular  cylinder  (without 
recourse  to  rediscretization).  They  used  Pohlhausen's  boundary-layer  ap- 
proximation method  to  determine  the  position  of  the  separation  point  and 
the  vorticity  to  be  introduced  at  each  time  interval.   Their  analysis  pre- 
dicted a  separation  angle  of  9  =64  degrees,  which  is  considerably  smaller 
than  the  commonly  accepted  experimental  value  of  8  =81  degrees. 

The  question  concerning  the  decay  of  vortices  has  remained  practically 
obscure.   Although  several  researchers  employed  methods  to  cancel  opposite- 
ly-signed vortices  whenever  they  came  within  a  prespecified  distance  of 
each  other,  no  systematic  investigation  of  the  effect  of  decay  of  vorticity 
was  conducted.  Often  recourse  was  made  to  Prandtl's  suggestion  that  indi- 
vidual vortices  in  a  vortex  street  retain  only  about  60  percent  of  the 
vorticity  generated  in  the  boundary  layer  during  a  period  required  to  shed 
a  single  vortex.   Some  investigators,  in  efforts  to  obtain  results  in 
partial  conformity  with  Prandtl's  suggestion,  automatically  applied  a  40 
percent  reduction  to  the  strength  of  each  nascent  vortex.   Seldom  was  men- 
tion made  of  the  experimental ly  observed  fact  that  the  rate  of  decay  of 
vortex  strengths  varies  with  distance  behind  the  cylinder. 

The  foregoing  represents  the  state-of-the-art  prior  to  the  beginning 
of  the  present  investigation.   A  relatively  more  detailed  account  of  the 
historical  evolution  of  the  DVM  is  found  in  Clements  and  Maul  I  [48]  and 
Fink  and  Soh  [46]. 
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III.   THEORETICAL  AND  PHENOMONOLOGICAL  FOUNDATIONS 
OF  THE  DISCRETE  VORTEX  MODEL 

A.  INTRODUCTION 

This  chapter  deals  with  the  characteristics  of  inviscid  and  viscous 
flows  about  a  circular  cylinder  from  the  forward  stagnation  point  to  the 
vortices  in  the  far  wake.  Mathematical  principles,  boundary  layer  con- 
cepts, and  relevant  hypotheses  are  introduced  in  preparation  for  the  dis- 
cussion of  the  details  of  the  discrete  vortex  model  developed  in  the 
present  study. 

B.  FLOW  KINEMATICS 

The  complex  potential  for  a  two-dimensional  inviscid  flow  is  given  by 


w(z)  =  (J)  +  \\p 

where  $  and  ij)  represent,  respectively,  the  velocity  potential  and  the 
stream  function.   The  velocity  components  which  satisfy  the  Cauchy-Riemman 
conditions  are  given  by 

u  =  -  ■si  =  —  -»j-  and  v  =  -  -s^  -  •**■         (5) 
dx     dy  3y   3x 

In  terms  of  the  complex  potential,  one  has 

dw(z)       ,  .  .,s 

—t- —  =  -  u  +  iv  (6) 

dz 

The  flow  about  a  circular  cylinder  is  generated  by  combining  a  uniform 

flow  with  a  doublet  at  the  center  of  the  cylinder  (see  e.g.,  Milne-Thomson 
[49]).   This  yields 

w(z)  =  w  =  -  U(z  +  c2/z)  (7) 
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in  which  -Uz  represents  the  uniform  flow  and  the  remainder  the  doublet. 
The  introduction  of  a  vortex  into  the  flow  field  exterior  to  the  cylinder 
requires,  by  virtue  of  the  circle  theorem,  the  introduction  of  an  image 
vortex  interior  to  the  circle.  Thus,  the  flow  about  a  circular  cylinder 
with  N  number  of  real  vortices  is  given  by 

2     .   N  2 

w  =  -U(z  +  — )  +  4-  V*  r  Un(z-z  )  -  &n(z  -  £-)}         (8) 
z     2tt  ^^   n      n  z 

i  n 

n=  I 

It  should  be  noted  that  a  real  vortex  located  at  z  -  z  has  an  imaqe  at 

n 

2  - 

z  =  c  /z  and  at  z  =  0.   The  image  at  the  center  of  the  circle,  which 
n 

gives  rise  to  terms  like  TJln(z),  must  be  excluded  due  to  the  condition 

that  the  vortices  shed  from  the  cylinder  leave  a  circulation  opposite  to 

their  own  on  the  cylinder.   Clearly,  the  absence  or  presence  of  an  image 

vortex  at  the  center  of  the  circle  does  not  alter  the  boundary  condition 

(9(f)/3n|      =  0)  since  the  streamlines  of  a  vortex  are  circles.  No  re- 

|z|=c 
striction  has  been  imposed  regarding  the  time  dependence  of  the  velocity 

U  other  than  that  it  be  uniform. 

Substituting  Eq.  (8)  into  Eq.  (6),  one  obtains  the  velocity  at  any 

point  z  as 

2       N 
_u+iv  =  dw  =  _u(|   c)+  i  y  {_J 1    } 

dz         2   2tt  L-i        n  z-z        2 
z        .       n      c 
n=l  z  -  -=- 

z 
n 

Equation  (9)  may  also  be  used  to  evaluate  the  convective  velocity  of  a  vor- 
tex at  z  =  z  provided  that  the  singular  term  T   /(z-z  )  is  omitted  by  virtue 
n  3  n     n  ' 

of  the  condition  that  a  vortex  cannot  impart  a  velocity  on  itself. 
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It  is  advantageous  to  work  with  dimension  less  parameters  by  introduc- 
ing the  following  change  of  variables: 

z"  =  z/c,  u'  +  iv'  *  (u+iv)/U,  I"  =  T/Uc  (10) 

w»(z)  =  w(z)/Uc      t1  =  Ut/c 

Substituting  into  Eqs.  (8)  and  (9),  letting  U  =  1.0,  c  =  1.0,  and  dropping 
the  primes  for  sake  of  simplicity,  one  has 

N 

wCz)  -  -(z+i)  +  4-   y  r  Un(z-z  )    -  4n(z-l/z  )}  (II) 

z    2tt  Z—/   n       n  n 

n=l 
and 


N 

-u+iv  =   -L  -  i  +  ±  y   V  {-J L_} 

z2       27T  Z,  n  z-zn      /; 

n=l  n 


(  12) 


where  u,  v,  r  and  z  are  now  the  normalized  variables  and  w  the  normalized 
complex  potential  function.   Hereafter,  all  measures  of  distance,  time, 
velocity  and  circulation  are  given  in  terms  of  the  above  normalized 
variables.  Any  exceptions  are  clearly  specified. 

Equation  (12)  is  used  extensively  in  convecting  the  vortices  and  in 
calculating  the  lift  and  drag  forces.   In  fact,  about  90  percent  of  the 
computation  time  of  the  present  model  is  consumed  in  the  evaluation  of 
complex  velocities.  This  points  out  the  necessity  of  keeping  the  number 
of  discrete  vortices  as  small  as  possible  while  maintaining  a  sufficiently 
small  time  interval  to  insure  that  the  characteristics  of  flow  are  faith- 
fully reproduced.   The  selection  of  the  time  interval,  like  the  selection 
of  the  mesh  size,  is  one  of  the  subjective  items  of  the  method  and 
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requires  a  careful  parametric  analysis.  This  is  addressed  to  in  a  later 
section. 

In  general,  the  calculation  of  velocities  and  the  convection  of  vor- 
tices are  straightforward  and  noniterati ve.   Under  certain  circumstances 
vortex-vortex  proximity  between  adjacent  spirals  and  the  vortex-boundary 
proximity  may  require  special  attention.  The  resolution  of  such  secondary 
issues  is  deferred  to  another  section  since  they  are  not  of  any  major  con- 
sequence either  theoretically  or  practically. 

C.  FORMULATION  OF  THE  RESISTANCE  EQUATIONS 

The  forces  acting  on  a  stationary  circular  cylinder  of  unit  length  by 

a  time-dependent  flow  containing  N  vortices  of  strengths  V      and  locations 

z  can  be  calculated  from  the  normalized  form  of  the  generalized  Blasius 
n  3 

theorem  (Milne-Thomson  f49])  to  give 


cd+  !cl4  /(dT)2dz+  ;lt  /wdz  (l3) 

where  C_  and  C  represent  the  drag  and  lift  coefficients  respectively. 

Inserting  Eq.  (II)  into  Eq.  (13)  and  evaluating  the  residues,  Sarpkaya 

L~50n  has  shown  that 

N 

Cn  +  iC,  =  V*  T  {(iu  -v  )  -  (iu  '-v  ')}  (14) 

D     L    Lmd        n     n   n       n    n 

n=l 

in  which  u  and  v  represent  the  velocity  components  of  the  n-th  real  vor- 
tex and  u   and  v   ,  the  velocity  components  of  the  image  of  the  n-th  vortex. 
Equation  (14)  may  be  written  in  a  more  compact  and  computationally 

more  suitable  form  as 

N 

C  +  iC.  =i  V  r  (z  -  z' )  (15) 

U     L     i—d        n   n     n 

n=l 
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where  the  dot  denotes  derivative  with  respect  to  time  such  that 


z  =  q  =  u  +  i  v  (1 6a) 

n   ^n    n     n 

and 

•  i     i     i   .   i  , . -.  , 

z=q=u+iv  (1 6b) 

n     n     n     n 

The  complex  velocities  q  and  q   may  easily  be  shown  to  be  related  by 


or  by 


q 

q   = *■  (17a) 

^n        2 

z 
n 


-u  +  i  v 
i  .  .   i     n     n  ,  .  _,  , 

u   +  iv   =  * (17b) 

n      n       -2 
z 
n 


The  normalized  force  given  by  Eqs.  (14)  or  (15)  can  also  be  evaluated  in- 
directly through  the  use  of  surface  pressures  and  the  unsteady  form  of 
Bernoulli's  equation.   In  either  case  the  results  are  identical  and  the 
use  of  one  method  or  the  other  is  simply  one  of  personal  choice  and  pro- 
gramming efficiency. 

D.   GENERATION  AMD  INTRODUCTION  OF  VORTICITY 

The  vorticity  per  unit  length  of  a  thin  shear  layer  separating  two 
streams  with  velocities  U.  and  IL  is  equal  to  the  difference  of  the  velo- 
cities U.-IL.  This  vorticity  is  transported  with  the  average  velocity  of 
the  two  streams,  i.e.,  (U,  +  l_L)/2.   This  is  a  continuous  process  and  the 
amount  of  vorticity  transported  during  a  time  At  becomes 

AT  =  (U|-U2)(U[+U2)   At/2  =  (Uj  2-U22)/2  .At         (18) 
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The  discrete  vortex  model  represents  the  continuous  generation  and  flux 
of  vorticity  into  the  outer  flow  by  introducing  point  vortices  at  discrete 
time  intervals.  The  strength  and  position  of  introduction  of  each  vortex 
must  depend  on  the  interaction  of  the  unsteady  vorticity  distribution  in 
the  wake  and  on  the  appropriate  boundary- layer  considerations  necessary  to 
determine  where  the  separation  occurs.  Those  applications  of  the  DVM 
which  consider  the  separation  points  fixed  ignore  not  only  the  said  inter- 
action but  also  artificially  constrain  the  evolution  of  the  flow. 

In  the  present  analysis  the  instantaneous  position  of  the  separation 
points  is  determined  through  the  use  of  an  approximate  boundary  layer 
analysis.  At  a  given  instant  the  velocity  distribution  on  the  forebody 
of  the  cylinder  is  calculated  and  the  corresponding  pressure  distribution 
is  assumed  to  be.  impressed  on  the  boundary  layer.   This  fundamental  feature 
of  the  boundary  layer  theory  together  with  the  integral  momentum  equation 
of  von  Karman  permits  the  use  of  a  number  of  separation  prediction  methods. 
Among  the  most  notable  of  such  methods  are  those  given  by  Pohlhausen  C5l3, 
Stratford  C52],  Thwaites  C53],  Curie  and  Skan  [54],  Timman  [55],  and 
Loitsianski  [56].  All  of  these  methods  exhibit  varying  degrees  of  sensi- 
tivity to  the  external  velocity  distribution  and  may  not  all  yield  identical 
or  uniformly  consistent  results. 

It  is  important  to  note  that  the  combination  of  potential  and  viscous 
flow  methods  cannot  establish  a  one-to-one  relationship  between  the  Reynolds 
number  and  the  flow  characteristics.   It  is  only  the  kind  of  separation 
criteria  used  (laminar  or  turbulent)  that  determines  whether  the  flow  is 
subcritical,  critical,  or  turbulent.   The  fact  that  the  use  of  the  DVM  be- 
comes increasingly  more  accurate  with  decreasing  shear  layer  thickness 
makes  one  realize  that  the  results  obtained  with  a  laminar  separation 
criteria  correspond  to  moderate  to  high  subcritical  Reynolds  numbers. 


The  integra I -momentum  methods  such  as  that  of  Pohlhausen  require  that 
the  velocity  outside  the  boundary  layer  be  given  and  that  the  separation 
points  be  either  fixed  or  fluctuate  with  small  amplitudes  and  velocities 
about  a  mean  position.   In  other  words,  the  changes  in  the  outer  flow  and 
the  separation  points  should  be  such  that  the  flow  may  be  treated  as 
quasi-steady. 

The  outer  velocity  distribution  may  be  calculated  either  numerically 
through  the  use  of  an  interactive  scheme  between  the  wake  and  the  boundary 
layer  or  through  the  use  of  experimentally  determined  pressure  distribution, 
Hiemenz  C57~J  used  the  measured  pressure  distribution  about  a  circular 
cylinder  and  approximated  the  velocity  distribution  at  the  edge  of  the 
boundary  layer  by 

UC5)  =  1.8348?  -  0.27422C3  -  0.04782£5  (19) 

where  E,   represents  the  arc  length  measured  along  the  cylinder  from  the 
stagnation  point.   Hiemenz's  solution  indicates  that  separation  occurs  at 
an  angle  8  =  82  degrees,  measured  from  the  stagnation  point.   This  cele- 
brated result  has  often  been  taken  as  proof  that  in  subcritical  flows 
separation  occurs  at  9  =82  degrees.   The  fact  has  not  been  realized  that 
the  entire  velocity  distribution  and  in  particular  that  portion  between 
the  maximum  velocity  and  the  separation  point  (i.e.,  from  8  =  70  to  8  - 
80  degrees)  dramatically  affect  the  position  of  the  separation  point.   In 
fact,  Meksyn  C58~J  has  shown  that  the  higher  order  derivatives  of  Hiemenz's 
velocity  profile  do  not  satisfy  certain  fundamental  criteria  regarding 
the  separation  point.  At  present  there  is  no  completely  satisfactory  velo- 
city distribution  for  the  forebody  of  a  circular  cylinder  either  in  sub- 
critical  or  in  the  supercritical  regimes.   Additional  discussion  of  this 
subject  is  found  in  Meksyn  [5"8]. 
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Numerous  experiments  have  yielded  laminar  separation  points  ranging 
from  9  =  70  degrees  to  9  =90  degrees.  The  commonly  accepted  value  is 
about  80  degrees.  Fage  and  Johansen's  D  0]  experiments,  and  others  since 
then,  have  shown  that  the  separation  point  is  in  reality  a  separation  zone 
which  extends  over  a  finite  region  and  that  spanwise  and  chordwise  fluctua- 
tion of  this  region  do  not  permit  a  precise  determination  of  the  mean 
position,  amplitude,  and  period  of  oscillation  of  the  separation  point. 

The  purpose  of  the  foregoing  brief  discussion  is  to  point  out  both  the 
theoretical  and  experimental  difficulties  encountered  in  the  determination 
of  the  separation  point  and  the  inadequacies  of  the  appropriate  laminar 
boundary  layer  methods.   The  prediction  of  the  separation  of  turbulent 
boundary  layers  is  even  more  uncertain  (Stratford  C52J,  Takada  [59])  and 
are  not  discussed  herein.   Instead,  attention  is  focused  on  Pohlhausen's 
method  for  time-dependent  boundary  layers.  The  former  is  more  appropriate 
to  the  flow  conditions  where  the  motion  of  the  separation  point  and  the 
velocity  at  the  edge  of  the  boundary  layer  may  be  treated  as  quasi-steady. 
Schuh's  C60]  method  is  more  appropriate  for  the  initial  instants  of  motion 
where  the  separation  points  move  rapidly  from  9  =180  degrees  to  about 
9  =  110  degrees.   The  motion  of  the  separation  point  commences  after  a 
relative  fluid  displacement  of  (Schlichting  L6I]) 

+  =  2d  hm  =  0-551  (20) 

The  reason  for  this  is  that  in  impulsively  started  flows  about  rounded 
bluff  bodies  the  separation  does  not  begin  until  the  fluid  moves  some  dis- 
tance past  the  body.   The  only  exception  to  that  is  the  flow  about  bodies 
with  sharp  corners  in  which  case  the  separation  begins  at  t  =  0.0. 

For  the  case  of  impulsively  started  flow  about  a  cylinder,  the  time 
in  the  boundary  layer  calculations  is  measured  from  the  instant  of  the 
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inception  of  separation.  Thus,  a  relative  time  of  t  =  0.351  must  be  added 
to  that  used  in  the  boundary  layer  calculations  in  order  to  match  the 
starting  times. 

Pohlhausen's  method  is  well  described  in  the  literature  (Schlichting 
[61])  and  only  the  highlights  of  the  method  are  presented  herein.  Assuming 
a  fourth  order  polynomial  for  the  velocity  distribution  within  the  boundary 
layer,  Pohlhausen  has  shown  that  von  Karman's  integral  momentum  equation 
may  be  reduced  to 


dZ   F(K)      .     „     dU(?)  In,  , 

and     K  =  Z   _,-  (21  ) 


in  which 


dC   U(T)  d£ 


K=  (3T5-  -  94T  A  "  9012  A')2a  (22) 

with 

A=A2dU(£!  (23) 

v   d£ 

Other  variables  appearing  in  Eqs.  (21),  (22)  and  (23)  are 

U(£):   the  normalized  velocity  on  the  cylinder 

£:   the  arc  length  measured  along  the  cylinder  from  the 
stagnation  point 

5:   the  boundary  layer  thickness  normalized  by  the  cylinder 
radi  us 

V:   viscosity  normalized  by  the  ambient  velocity  and  the 
cylinder  radius  (sometimes  referred  to  as  the  unit 
Reynolds  number). 

The  values  of  K  and  F(K)  are  tabulated  in  terms  of  A  (Schlichting  [61]). 

Similar  tables  are  provided  for  other  boundary  layer  approximations  cited 

previously. 

integration  of  the  coupled  Eqs.  (21)  begins  at  the  front  stagnation 

point  and  proceeds  along  both  sides  of  the  cylinder  independently  until 
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the  separation  criteria  (A  -   -12)  is  satisfied  for  each  layer.  This  re- 
sults in  two  distinct  separation  points  which  in  general  are  not  at 
identical  angular  distances  from  the  mobile  stagnation  point. 

As  noted  earlier,  the  foregoing  method  is  applicable  only  to  the  rela- 
tively established  flow  where  the  velocity  at  the  edge  of  the  boundary 
layer  and  the  separation  point  do  not  significantly  vary  with  time.   During 
the  early  stages  of  motion,  however,  the  time  dependence  of  the  said  para- 
meters is  very  strong  and  must  be  treated  through  the  use  of  the  time 
dependent  integral  boundary-layer  equations.  Schuh  C60]  has  extended 
Pohlhausen's  method  for  the  general  calculation  of  the  unsteady  boundary 
layer.   The  details  of  his  method  are  not  presented  here.   The  results 
obtained  by  Schuh  for  an  impulsively  started  flow  about  a  circular  cylin- 
der are  shown  in  Fig.  8.   These  results,  based  on  the  classical  potential 
flow  velocity  distribution  about  a  circle,  show  that  the  separation  point 
rapidly  approaches  a  value  of  about  0  =  109  degrees.  Schuh's  method  is 
not  capable  of  incorporating  the  growth  of  the  symmetric  wake  and  the 
subsequent  motion  of  the  separation  points  and  thus  is  not  valid  at  large 
times.   Thoman  and  Szewczyk  C20]  employed  the  finite  difference  form  of 
the  full  Navier-Stokes  equations  and  obtained  results  identical  to  those 
shown  in  Fig.  8  for  the  initial  instants  of  flow.   They  have  further  shown 
that  the  separation  point  once  havi ng  reached  an  angular  position  of  about 
MO  degrees  moves  gradually  to  the  final  value  of  the  separation  angle 
(about  80  degrees).   Partly  to  confirm  this  conclusion  and  partly  to  assess 
the  significance  of  time  dependence  on  the  instantaneous  position  of  the 
separation  points,  two  sets  of  calculations  are  performed.   Both  use  the 
time-dependent  velocity  distribution  available  through  the  model  and  both 
begin  with  an  initial  separation  point  of  0  =  109  degrees.   The  first  cal- 
culation uses  Pohlhausen's  steady-state  separation  criteria  with  the  in- 
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Fig.  8.   Separation  angle  versus  time  [603. 
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stantaneous  velocity  distribution.  The  second  calculation  uses  the  unsteady 
integra I -momentum  equations  given  by  Schuh  C6CG,  including  those  terms 
associated  with  the  unsteady  velocity  distribution.  The  results  obtained 
from  both  sets  of  calculations  are  shown  in  Fig.  9.   Evidently  the  two  cal- 
culations yield  almost  identical  results,  pointing  out  the  fact  that  the 
separation  point  moves  rather  slowly  and  that  the  contribution  of  the  time 
dependent  terms  in  the  unsteady  form  of  the  equations  are  negligible.   It 
is  because  of  this  reason  that  all  calculations  are  based  on  Pohlhausen's 
separation  criteria  assuming  the  flow  to  be  quasi-steady.  This  assumption 
improves  as  one  approaches  steady  state.   Figure  9  indicates  that  the 
angular  speed  of  the  separation  point  is  about  10  degrees  per  second.  The 
same  speed  about  the  mean  position  of  the  separation  point  is  about  2  de- 
grees per  second.   (As  is  shown  later,  the  separation  point  oscillates 
about  a  mean  position  of  9  =77  degrees  with  an  amplitude  of  3  degrees 
and  a  period  of  about  10.)  Based  on  the  foregoing,  an  amplitude  of  oscil- 
lation as  much  as  15  degrees  can  be  performed  without  the  inclusion  of 
time-dependent  terms. 

In  the  present  model  the  calculations  begin  at  the  time  of  t  =  1.0 
when  the  separation  point  is  at  8  =  109  degrees.   It  should  be  noted 
that  a  straightforward  application  of  Pohlhausen's  method  with  a  potential 
flow  velocity  distribution  (U  =  2  sin  8)  yields  a  separation  angle  of 
8  -    109.5  degrees. 

Once  the  position  of  the  separation  points  is  determined,  the  rate  of 
generation  of  vorticity  is  easily  determined  by  writing 


dr    C     3u 

dt-J  I7U 


dy  (24) 
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Fig.  9.   Comparison  of  steady  and  unsteady 
separation  point  calculations. 
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where  5  is  the  boundary  layer  thickness;  u,  the  velocity  within  the  bound- 
ary layer;  and  y,  the  distance  along  the  outward  normal.  The  contribution 
of  dv/dx    is  excluded  from  the  vorticity  expression  in  conformity  with  the 
approximations  of  the  boundary  layer  theory. 

The  integration  of  Eq.  (24)  to  obtain  the  circulation  flux  at  a  separa- 
tion point  yields 

£«V  (25) 

dt    2  ^D1 

where  U   is  the  tangential  velocity  at  the  predicted  separation  point. 
Ordinarily  the  rate  of  change  of  circulation  should  be  written  as 

2   2 
d"  _  I 2  ,9,. 

dt 2 (26) 

where  U,  and  IL  represent  respectively  the  velocities  at  the  outer  and 
inner  edges  of  the  separating  shear  layer.   Extensive  experiments  with 
various  bluff  bodies  by  Fage  and  Johansen  0  Oj  have  shown  that  IL  is  about 
5  percent  of  U,  and  that  -~  is  correctly  represented  by  U  12.      This  is  a 
universally  accepted  experimental  and  theoretical  fact. 

In  numerical  calculations  -rr  is  replaced  by  AIVAt  and  the  strength  of 
the  nascent  vortex  is  calculated  from 

U  2At 

r   =  -2* —  (27) 

nv     2 

The  placement  of  nascent  vortices  into  the  flow  to  adequately  represent 
the  separating  shear  layers  has  been  a  source  of  difficulty.   This  i-s  pri- 
marily because  the  discretization  of  a  continuous  process  near  a  singular 
region  (separation  point)  is  an  exceedingly  complex  problem.   In  the  present 
analysis,  no  new  procedures  are  needed  for  the  determination  of  the  strength 
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and  position  of  the  nascent  vortex  because  the  separation  point  is  obtained 
from  the  boundary  layer  analysis.  By  virtue  of  its  definition,  the 

separation  point  must  be  the  point  at  which  both  the  velocity  and  the 

2 
shear  stress  on  the  cy I inder  are  zero  .  Thus,  the  nascent  vortex  of  known 

strength  must  be  placed  on  a  radial  line  passing  through  the  separation 

point  at  a  distance  m,  i.e.,  at  z   =  (l+e)  exp  C I C tt— 9  )],  to  render  the 

tangential  velocity  on  the  cylinder  at  the  separation  point  equal  to  zero. 

This  yields  a  unique  value  for  e  as 

i  +  |r  |/2itu 

1  n  v '     s 
e  =  III 2.  -  |  (28) 

I  -  IT   |/2itU 
1  nv1    s 

The  use  of  the  boundary  layer  theory  and  the  separation  condition  uniquely 
specifies  the  strength  and  position  of  the  nascent  vortices  and  allows  the 
separation  points  to  move  along  the  cylinder  as  required  by  the  interaction 
between  the  wake  and  the  boundary  layer.   It  is  of  importance  to  note  that 
the  no-slip  condition  is  satisfied  only  at  discrete  time  steps.  Efforts 
to  satisfy  the  said  condition  in  a  more  continuous  manner  are  counterproduc- 
tive because  the  decrease  in  At  calls  for  unacceptably  large  computation 
times.  One  of  the  essential  features  of  the  model  is  the  bracketing  of  a 
suitable  interval  for  At,  which  provides  an  adequate  scale  of  discretiza- 
tion without  excessive  computation  time. 

The  method  described  above  properly  accounts  for  the  vorticity  generated 
on  the  forebody  of  the  cylinder.  However,  the  flow  in  the  wake  of  the 
cylinder  develops  a  boundary  layer  on  the  afterbody  which  generates 


"2 
A  lengthy  discussion  may  be  made  regarding  the  difference  between  the 

definitions  of  steady  and  time  dependent  separation  criteria.   However,  a 

numerical  analysis  such  as  the  DVM  does  not  warrant  the  consideration  of 

the  distinction  between  the  results  stemming  from  the  use  of  different 

separation  criteria. 
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countervorticity.  Procedures  by  which  this  relatively  small  amount  of 
vorticity  generation  may  be  accounted  for  are  far  from  simple.  The  velo- 
city distribution  on  the  afterbody  has  strong  spatial  and  temporal  gra- 
dients due  to  the  immediate  proximity  of  the  fluctuating  wake.  These 
gradients  render  useless  the  separation  prediction  methods  described 
earl ier. 

Davis  C62]  attempted  to  account  for  countervorticity  effects  by  using 
four  separation  points  (two  for  the  forebody  and  two  for  the  afterbody 
boundary  layers)  but  his  method  required  the  assumption  of  fixed  separation 
points.  Deffenbaugh  and  Marshall  C47H  approximated  the  rear  shear  layer 
effect  by  assuming  that  the  boundary  layers  remain  attached  over  the  same 
fraction  of  a  region  of  adverse  pressure  gradient  as  that  of  the  forward 
boundary  layers.   However,  in  view  of  the  complex  velocity  distribution 
on  the  afterbody,  it  is  incorrect  to  apply  the  same  boundary-layer  con- 
siderations as  were  used  on  the  forebody.   It  is  safe  to  say  that  the 
numerical  perturbations  introduced  into  the  model  in  attempting  to  account 
for  the  secondary  effects  of  countervorticity  may  cause  more  error  than 
would  the  neglect  of  the  effect  altogether.  This  consideration  plus  the 
introduction  of  the  circulation  reduction  mechanism  (yet  to  be  described) 
obviate  the  need  for  the  separate  accounting  of  the  effect  of  the  rear 
shear  layers. 

E.  CONVECTION  AND  REDI SCRETI ZATION  MECHANISM 

The  vortices  emanating  from  a  separation  point  may  be  convected  in  a 

number  of  ways.  Three  most  promising  schemes  are  those  given  by  the 

following  expressions: 

zCt+At)  =  z(t-At)  +  2q(t)  At  (29a) 

zCt+At)  =  z(t)  +  iC3q(t)-q(t-At)]  At  (29b) 

z(t+At)  =  z(t)  +  q(t)  At  (29c) 
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where  q(t)  is  the  complex  velocity  of  the  vortex  evaluated  at  the  positions 
and  times  indicated.  These  schemes  yield  more  or  less  the  same  accuracy 
and  the  use  of  one  in  lieu  of  the  others  is  primarily  governed  by  the  sim- 
plicity of  the  method  and  computation  time  considerations.   In  the  present 
analysis  the  simple  Eulerian  scheme  given  by  Eq.  (29c)  is  used  because  of 
its  simplicity  and  because  the  use  of  relatively  small  time  intervals  (At  = 
0.125)  make  the  method  sufficiently  accurate.  Clearly,  the  results  obtained 
with  the  three  schemes  cited  above,  or  with  others,  cannot  be  compared 
without  taking  into  consideration  the  effect  of  the  time  intervals  used 
and  the  discretization  process  (if  one  is  employed).  All  such  convection 
schemes  must  assure,  within  the  overall  limitation  of  the  model,  that  the 
vortex  paths  in  regions  of  very  small  radius  of  curvature  do  not  signifi- 
cantly deviate  from  those  calculated  with  much  smaller  time  intervals. 

It  has  been  explained  in  connection  with  the  discussion  of  the  histori- 
cal development  of  the  DVM  that  all  investigators,  with  the  exception  of 
Fink  and  Son  C46H,  let  the  vortices  move  in  the  form  of  vortex  clouds  and 
took  special  measures  to  eliminate  the  difficulties  arising  from  vortex- 
vortex  and  vortex-cylinder  proximity  effects.  Not  to  be  entirely  critical, 
the  use  of  vortex  clouds  can  yield  force-transfer  coefficients  nearly 
identical  to  those  obtained  with  the  use  of  redi screti zed  sheets.  However, 
one  drawback  of  the  vortex  clouds  is  that  it  does  not  allow  the  formation 
of  spiralling  vortex  sheets  as  often  observed  experimentally.  There  are 
other  features  of  the  red i screti zation  which  will  become  clear  later. 
The  vortex  clouds  stem  partly  from  the  accumulation  of  the  logarithmic 
errors  and  do  not  in  any  way  represent  a  turbulent  vortex.   In  the  early 
stages  of  the  present  investigation  the  vortex  clouds  were  used  to  obtain 
some  idea  about  the  gross  features  of  the  asymmetric  wake.   Figure  10  is 
a  sample  computer  plot  of  vortex  clusters.   It  is  not  possible  to  identify 
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Fig.  10.   Flow  configuration  without  red i scret i zat ion . 


Fig.  II.   Flow  configuration  with  red i scret i zat ion . 
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a  spiralling  sheet  in  Fig.  10,  but  the  shedding  of  the  vortices  and  the 
thinning  of  the  connecting,  vortex  sheets  between  two  adjacent  vortex 
clusters  are  clearly  observable. 

The  fact  that  the  vortices  must  represent  each  segment  of  the  shear 
layer  as  accurately  as  possible  and  thus  avoid  the  buildup  of  logarithmic 
errors  in  their  convection  gave  rise  to  the  method  of  rediscreti zation  or 
redistribution.  A  sample  flow  pattern  obtained  with  rediscreti zation  is 
shown  in  Fig.  II  where  the  point  vortices  are  connected  by  line  segments 
to  denote  the  position  of  the  sheet.  The  mathematical  foundations  of  the 
concept  have  already  been  introduced.   In  the  following,  the  details  of 
the  mechanism  where  the  rediscretization  is  affected  are  described. 

Consider  a  segment  of  the  vortex  sheet  at  an  instant  t  and  assume  that 
the  vortices  are  connected  by  straight  line  segments  (see  Fig.  12).   The 
distance  between  any  two  vortices  is  given  by 


Z4 


Fig.  12.  Vortex  sheet  before  rediscretization, 


6s  =  z  -  z  (30) 

n   '  n    n- I 

The  length  of  the  sheer  from  the  position  z.  (core  of  the  spiral)  to  the 
n-th  vortex  is  given  by 
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n 


n 


•«  ■  E  ^  ■  23 


z  .  -  z  . 

J    ^— '   '  J    J" 

j=2       J-2 


(3D 


The  total  length  of  the  sheet  for  N  number  of  vortices  becomes 


N 


S  5s 


N    A-rf   "J 

J-2 


(32) 


The  vorticity  per  unit  length  of  the  sheet  at  each  z   is  given  by 


Yn  =  (6s  ^,-6s   ~)77 
n+l    n-l 


YM  =  r./Ss.. 

N    N   N 


N  <  n  <  I 


(33) 


Y|  =  r,/6s2 


The  circulation  per  unit  length,  y   ,  may  be  tabulated  as  a  function  of  s 
through  the  use  of  Eqs.  (31)  and  (33).  This  tabulation  is  then  used  as  a 
piecewise  continuous  function  denoted  by  y  -  y(.s) .      The  function  y(s)  uses 
a  linear  interpolation  of  the  tabulated  values,  i.e. 


Y  —  Y 

r   i       , i  \  n  n-l 

y(s)  =  y  .  +(s-s   ,)  

'      'n-l      n-l   s  -s   . 

n  n-l 


s   ,  <  s  <  s 

n-l  —   —  n 


(34) 


YCs)  represents  circulation  per  un i f  length  and  is  an  invariant  at  a  given 
time  of  the  red iscret i zation  process. 

The  next  step  in  the  analysis  is  the  replacement  of  the  existing  un- 
evenly  spaced  N  vortices  with  N  evenly  spaced  new  vortices.   In  general  N 
may  be  different  from  N.   Replacement  of  N  vortices  is  accomp I  i shed  by 
starting  at  the  core  of  the  spiral,  z.,  and  placing  vortices  at  intervals 
of  5s  measured  along  the  straightline  segments  connecting  the  original 


65 


vortex  positions  (see  Fig.  13).   It  should  be  noted  that  when  the  method  is 
applied  to  the  example  depicted  in  Fig.  13,  z„  and  z,  are  positioned  such 


that  z. 


=  6s  and  | z,  -  zJ  +  |z2 


z-|  =  6s.   In  other  words,  dis- 


"2    I 
tances  are  measured  along  the  straightline  segments  of  the  existing  sheet. 


Fig.  13.   Vortex  sheet  after  red i sc re t i za t i on  . 

The  strengths  of  the  new  vortices  are  calculated  through  the  use  of  the 
distribution  function  y(s)    such  that  the  circulation  per  unit  length  at  a 
new  position,  z  ,  is  the  same  as  for  the  original  sheet  at  the  said  posi- 
tion, i.e., 


F.  =  y(.o)    6s, 


y(s*) 


6s  ^,+6s 
n+l   n-l 


<  n  <  N 


(35) 


where 


Y(sN)  6s 


s*  =  (n-l)  6s 

The  final  step  in  the  basic  procedure  is  to  insure  that  there  has  been  no 
net  gain  or  loss  of  circulation  during  the  red i scretization  process.   The 
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circulation  deficit  is  computed  as 

N        N 

Ar  =  V  r  -  V  r  (36) 

'    '        n    /  '        n 
n=l      n»l 

and  an  equal  portion  of  the  difference  is  added  algebraically  to  each  vor- 
tex, i.e., 

AT 

r  =  r  +  —  (37) 

n    n 

N 

Two  parameters  yet  to  be  specified  are  6s  and  N.  The  present  model 
does  not  change  the  number  of  vortices  during  red i scretization  and  thus 
N  =  N.  6s  and  N  are  related  by 


6s  =  s  /(N  -  I  )  (38) 


It  is  important  to  note  that  N  is  the  total  number  of  vortices  in  a  detached 
sheet  (The  definition  of  a  detached  sheet  is  introduced  later.)  and  the 
redi screti zat ion  process  does  not  change  the  position  of  either  the  first 
or  the  N-th  vortex.   For  the  vortex  sheets  attached  to  a  nascent  vortex, 
however,  N  represents  the  number  of  vortices  excluding  the  nascent  vortex. 
The  reason  for  this  is  that  the  conditions  (strength  and  location  of  the 
nascent  vortex)  prevailing  at  or  near  the  separation  point  should  not  be 
altered  by  any  process  once  the  nascent  vortex  has  been  placed  in  the  flow 
in  accordance  with  the  no-slip  condition.  Thus,  for  an  attached  sheet  the 
red i screti zat ion  process  stops  at  the  vortex  immediately  adjacent  to  the 
nascent  vortex. 

It  is  appropriate  at  this  time  to  give  some  idea  about  the  actual 
magnitude  of  some  of  the  parameters  mentioned  above: 
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(i)    the  strength  of  a  nascent  vortex  -   0.05  to  0.15 
(ii)   N  =  3  to  80  depending  on  the  length  of  the  sheet 
Ciii)  6s  =  0.1 

The  above  numbers  are  representative  of  those  encountered  for  a  stationary 

cylinder  with  a  time  step,  At  =  0.125. 

F.   VORTEX  SHEDDING  AND  CIRCULATION  REDUCTION 
I .  Vortex  Shedding 

One  of  the  least  understood  of  the  wake  formation  processes  is  vor- 
tex shedding  where  the  sheet  connecting  the  spiralling  vortex  to  the  sepa- 
ration point  is  cut.   That  portion  of  the  sheet  which  remains  attached  to 
the  cylinder  begins  to  roll  up  into  a  new  vortex  spiral  while  the  detached 
portion  moves  downstream  to  become  part  of  the  vortex  street.   The  shedding 
process  is  accompanied  by  a  number  of  interactive  phenomena  such  as  the 
motion  of  the  shed  vortices,  adjustment  of  the  longitudinal  and  lateral 
spacing  of  the  vortex  cores,  and  the  reduction  of  circulation  through  the 
exchange  of  oppositely  signed  vorticity.   These  phenomena  are  so  inter- 
locked that  attempts  to  explain  the  reasons  for  the  occurrence  of  one  be- 
come in  turn  questions  to  be  explained.   Evidently  one  possible  way  to 
resolve  this  dilemma  is  to  have  an  exact  solution  of  the  problem.   Numeri- 
cal models  such  as  the  finite  difference  technique,  having  failed  to 
elucidate  the  interactive  mechanisms,  left  researchers  with  the  need  to 
develop  techniques  which  could  be  used  in  more  phenomonologica I  models 
such  as  the  DVM.   It  is  absolutely  essential  that  one  begin  the  analysis 
with  a  model  and  then  examine  its  validity  through  its  predictions. 

Inspired  by  flow  visualization  and  wake  measurements,  Gerrard  £.632 
proposed  that  the  growing  vortex  spiral  draws  the  opposite  shear  layer 
across  the  wake  when  it  becomes  sufficiently  strong.   The  approach  of 
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oppositely  signed  vorticity  cuts  off  further  supply  of  circulation  to  the 
vortex.  The  approximate  time  of  the  cutting  of  the  sheet  is  then  considered 
as  the  beginning  of  the  shedding  process.  This  mechanism  is  sketched  in 
Fig.  14. 

The  wake  comprised  of  such  shed  vortices  continuously  adjusts  it- 
self in  both  lateral  and  transverse  spacing.   In  general  if  is  agreed  that 
the  wake  may  be  divided  into  three  parts.  The  first  is  the  so  called  for- 
mation region  (£,  -   5);  the  second,  the  stable  region  where  the  vortices 

exhibit  the  character isttcs  of  a  fairly  uniform  vortex  street  (5  <  £   < 

'  —  sr  — 

12);  and  the  third,  the  unstable  region  extending  beyond  x  larger  than 
about  12  [64]. 

The  lateral  spacing  of  the  vortex  cores  gradually  increases  within 
the  stable  region  and,  like  the  longitudinal  spacing,  depends  on  the 
Reynolds  number. 

Even  though  considerable  experiments  have  been  conducted  (particu- 
larly at  low  Reynolds  numbers)  to  delineate  the  basic  features  of  the  wake, 
very  few  attempts  have  been  made  to  measure  directly  the  strength  of  the 
vortices  as  a  function  of  x  at  a  given  Reynolds  number.   Admittedly,  this 
is  an  extremely  difficult  undertaking.   Several  investigators  (Schaefer 
and  Eskinazi  [64],  Griffin  [65],  and  Davies  [66])  measured  the  longitudinal 
and  transverse  spacing  of  the  vortices  in  the  stable  region  and  calculated 
indirectly  the  strength  of  the  vortices  (assumed  to  be  identical)  through 
the  use  of  von  Karman's  vortex  street  model.   Consequently,  the  vortex 
strengths  are  approximate  and  highly  dependent  on  the  initial  inputs.   The 
most  important  conclusion  resulting  from  such  investigations  is  that  the 
strength  of  the  vortices  decreases  with  increasing  distance  downstream. 
The  question  concerning  how  much  of  the  circulation  generated  remains  in 
a  vortex  has  been  the  subject  of  discussion  since  the  time  of  Prandtl.  As 
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noted  earlier  Prandtl  suggested  that  the  net  circulation  remaining  in  a 
vortex  is  about  50  percent  of  that  generated.  Discrete  vortex  models  used 
by  Clements  C67]  and  Sarpkaya  C68H  give  a  reduction  of  about  15  percent. 
This  value  i s  s ign i f icant I y  smaller  than  the  estimates  of  40  percent  re- 
duction (Fage  and  Johansen  C ' 0~J ) ,  57  percent  (Roshko  [J  Q)  and  70  percent 
CBIoor  and  Gerrard  C69H).   It  should  not  be  inferred  from  the  foregoing 
that  the  vortices  once  having  acquired  a  net  circulation  of  about  40  per- 
cent continue  to  retain  that  circulation  forever.   The  understanding  of 
the  physics  of  circulation  reduction  is  of  major  importance  in  the  devel- 
opment of  any  phenomono logi ca I  model  and  is  discussed  in  greater  detail  in 
the  following  section. 

2.  Circulation  Reduction 

All  investigators  using  the  DVM  have  recognized  that  greater  atten- 
tion must  be  paid  to  the  circulation  reduction  if  the  DVM  is  to  produce 
results  comparable  to  those  obtained  experimentally.   In  the  early  stages 
of  the  evolution  of  the  present  model  the  need  for  a  circulation  reduction 
mechanism  was  quickly  reaffirmed  from  the  magnitudes  of  the  iift  and  drag 
forces. 

Before  undertaking  a  detailed  discussion  of  possible  circulation 
reduction  mechanisms,  it  is  advantageous  fo  discuss  the  results  obtained 
through  the  measurement  of  the  vortex  strengths.   Schmidt  and  Tilmann  [70], 
through  the  use  of  an  ingenious  ultrasonic  sound  beam  technique,  measured 
the  variations  of  the  circulations  of  vortices  without  inserting  a  probe. 
Their  measurements  were  made  at  various  Reynolds  numbers  by  changing  both 
the  velocity  of  flow  and  the  diameter  of  the  cylinder.   Figure  15  shows  T 
as  a  function  of  x  and  is  representative  of  the  results  obtained  by  Schmidt 
and  Tilmann.   In  their  experiments  the  Strouha I  number  was  about  0.21  and 
the  mean  velocity  at  the  point  of  separation  may  be  taken  as  U  -  I .5  U. 
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Fig.  I5.   Vortex  strength  versus  distance  [69,  70], 
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Fig.  16.   Percent  of  circulation  reduction  versus  distance  [69,  70]. 
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Thus,  the  total  circulation  generated  during  the  period  T  =  l/f  is  approxi- 
mated by 

r  =iu  2-°—=  10  7 

1    2  s  U  St   IU-/ 
o 

In  the  observed  formation  region,  I  <  x  <  6,  the  vortices  grow  rapidly  and 
at  the  same  time  lose  some  of  their  strength  due  to  various  circulation  re- 
duction mechanisms.  One  can  think  of  a  net  circulation  reduction  on i y  for 
vortices  which  are  in  the  stable  and  unstable  regions.   It  must  be  noted 
that  the  ultrasonic  technique  used  in  Schmidt  and  Tilmann's  experiments 
cannot  yield  accurate  vortex  strengths  in  the  formation  region  because  the 
path  of  the  sound  beam  crosses  over  two  oppositely  signed  vortices.   In 
the  stable  region,  however,  only  one  vortex  crosses  the  path  of  the  sound 
beam  at  a  given  time.   For  the  reasons  cited  above,  part  of  the  curve  ex- 
tending from  x  -  1 0  to  its  experimentally  reported  maximum  value  of  40  is 
used  to  calculate  the  circulation  reduction  by  dividing  the  circulation 
values  by  10.7.   The  result  is  shown  in  Fig.  16.   Also  shown  in  Fig.  16 
are  two  reduction  ratios  computed  similarly  from  Bloor  and  Gerrard's  E69] 
reported  vortex  strengths  at  x  values  of  12  and  20.   The  circulation  re- 
duction shown  in  Fig.  16  should  be  interpreted  only  as  a  reflection  of 
the  change  of  circulation  with  distance  rather  than  as  definitive  values 
since  such  values  depend  on  the  Reynolds  number,  three-dimensionality  ef- 
fects, and  the  methods  of  their  evaluation. 

One  or  more  physical  mechanisms  are  needed  to  explain  and  to  model 
the  observed  circulation  reduction.  Circulation  may  be  reduced  only  by 
oppositely  signed  circulation.   This  hypothesis  is  somewhat  mechanistic 
but  is  easily  adaptable  to  numerical  modelling.   There  is  no  question  that 
in  reality  the  entire  mechanism  may  be  related  to  the  action  of  viscosity. 
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Leaving  aside  momentarily  the  complex  relationship  between  the  action  of 
viscosity  and  circulation  reduction,  attention  is  concentrated  on  three 
sources  of  circulation  exchange.   It  has  been  noted  by  the  users  of  the 
DVM  that  the  circulation  may  be  reduced  by  the  oppositely-signed  circula- 
tion generated  on  the  afterbody  (countervorticity) .   In  accordance  with 
this  assumption,  the  discrete  vortices  which  come  closer  than  a  preassigned 
distance  to  the  cylinder  and  the  real  vortices  which  cross  the  cylinder 
boundary  are  annihilated.  Experiments  by  Fage  and  Johansen  DO]  have  shown 
that  only  about  a  10  percent  reduction  is  brought  about  by  the  entrainment 
of  fluid  bearing  circulation  of  opposite  sign  from  the  flow  aft  cf  the 
separation  points.   Numerical  experiments  to  be  discussed  in  another  sec- 
tion of  this  work,  as  well  as  those  done  by  others,  show  that  the  circula- 
tion generated  on  the  downstream  face  of  the  cylinder  can  account  for  no 
more  than  about  1/5-th  of  the  total  reduction  in  circulation.   This  fact 
calls  for  other  mechanisms  which  account  for  the  remainder  of  the  circula- 
tion loss.   The  second  mechanism  postulated  for  the  exchange  of  oppositely- 
signed  circulation  is  an  integral  part  of  the  shedding  mechanism.   As  pro- 
posed by  Gerrard  L"633,  when  the  growing  vortex  acquires  a  sufficiently 
large  strength,  it  draws  the  feeding  shear  layer  of  the  other  growing  vor- 
tex across  the  wake,  thereby  inhibiting  further  growth  of  the  vortex. 
The  tail  of  the  cut  sheet  is  then  absorbed  by  the  still  attached  vortex  on 
the  opposite  side  of  the  wake.  This  mechanism,  quite  distinct  from  the 
effect  of  the  boundary  layer  on  the  rear  of  the  cylinder,  further  reduces 
the  strength  of  a  given  vortex.   However,  numerical  experiments  again  show 
that  the  reduction  due  to  this  mechanism  alone  is  far  from  sufficient  to 
bring  the  circulation  down  to  the  observed  or  measured  values.   Consequent- 
ly, a  third  mechanism  must  be  postulated  with  which  not  only  the  vortices 
in  the  formation  region  but  also  those  further  downstream  experience 
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additional  circulation  reduction.  The  origins  of  the  third  reduction  mech- 
anism may  be  partly  in  the  ingestion  or  entrainment  of  fresh  fluid  into  the 
shear  layers,  partly  in  the  instabilities  in  the  spiralling  vortex  sheets, 
partly  in  slow  viscous  reduction,  and  partly  through  some  other  as  yet 
unknown  mechanism.  Regardless  of  the  true  origins  of  the  reduction  mech- 
anism, two  facts  remain  undisputed:   the  circulation  of  a  given  vortex  does 
decay  with  time  and  the  vorticity  diffuses  through  viscous  and  turbulent 
action;  and,  the  DVM  cannot  predict  even  the  most  gross  features  of  the 
real  flows  (for  example,  lift  and  drag  coefficients  and  circulation  distri- 
butions) without  the  use  of  a  circulation  reduction  mechanism. 

Lastly,  it  is  necessary  to  briefly  discuss  the  effect  of  diffusion 
of  vorticity  on  the  apparent  reduction  of  circulation.   Diffusion  means 
that  the  vorticity  is  spread  over  a  larger  area.   This  in  turn  means  that 
larger  numbers  of  suitably  distributed  point  vortices  for  each  segment  of 
the  vortex  sheet  rather  than  a  single  elemental  vortex  represent  the  corre- 
sponding distributed  vorticity.  Ordinarily,  the  point  vortices  (whether 
on  a  sheet  or   concentrated)  are  assumed  to  be  placed  at  the  center  of 
gravity  of  the  vorticity  field.   This  does  not  imply  that  the  single  vor- 
tex at  the  center  of  gravity  will  yield  the  same  velocity  at  an  arbitrary 
.point  on  the  cylinder  or  at  the  position  of  another  point  vortex  as  the 
distributed  vorticity.   As  far  as  the  cylinder  is  concerned,  the  said 


The  center  of  vorticity  of  M  point  vortices  is  computed  as 
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effect  becomes  relatively  small  as  the  distance  between  the  vortex  and  the 
cylinder  increases.  Consequently,  the  calculations  are  not  noticeably 
affected  by  the  diffusion  of  vorticity  in  the  far  field.  Since  diffusion 
is  a  matter  of  time,  the  vortices  in  the  immediate  vicinity  of  the  cylinder 
do  not  undergo  significant  diffusion.  The  foregoing  simply  points  out  the 
difficulty  of  accounting  for  all  direct  or  indirect  decay  mechanisms  and 
the  need  for  a  suitable  hypothesis  which  will  account  for  the  circulation 
decay  and  yet  be  compatible  with  the  basic  premises  of  the  DVM.  As  is 
discussed  in  greater  detail  later,  the  present  model  accounts  for  the  effect 
of  rear  shear  layer,  circulation  cancellation  brought  about  by  the  shedding 
mechanism,  and  finally  by  the  cross-street  circulation  exchange.  The  last 
process  is  accomplished  simply  by  transferring  a  certain  percentage  circula- 
tion between  oppositely  signed  vortices  at  each  time  step. 
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IV.   DETAILS  OF  THE  NUMERICAL  MODEL  AND  DISCUSSION  OF  RESULTS 

A.  SPECIFIC  DETAILS  OF  THE  MODEL' 

I .  Typical  Sequence  of  Calculations  at  an  Arbitrary  Time 

The  most  straightforward  way  to  describe  the  particulars  of  the 
model  is  to  examine  the  sequence  of  events  that  occur  during  a  typical  com- 
putational cycle.   Figure  17  is  a  flowchart  indicating  the  order  in  which 
the  model  executes  the  various  phases  of  the  analysis.   Each  of  the  events 
that  occur  during  every  computation  cycle  are  described  in  the  order  of 
their  occurrence.   The  asymmetry  introduction  and  vortex  shedding  processes 
occur  only  under  certain  flow  conditions  and  are  addressed  separately. 

At  the  beginning  of  each  computation  cycle  the  velocity  distribu- 
tion on  the  forebody  of  the  cylinder  is  calculated  at  one  degree  intervals 
from  the  most  recent  position  of  the  forward  stagnation  point.   The  posi- 
tions of  the  upper  and  lower  separation  points  are  then  determined  through 
the  use  of  Pohlhausen's  method  as  described  earlier.   This  allows  the 
strength  and  positions  of  the  two  nascent  vortices  to  be  determined  accord- 
ing to  the  no-slip  condition. 

The  next  step  is  the  convection  of  each  of  the  vortices  by  a  simple 
Eulerian  method  where  the  velocity  of  each  vortex  is  computed  at  its  pres- 
ent location.  The  velocities  of  the  nascent  vortices  are  computed  at  the 
point  of  their  introduction  prior  to  their  actual  placement  into  the  flow. 
This  insures  that  a  nascent  vortex  is  convected  with  the  velocity  prevail- 
ing at  the  time  and  place  of  its  creation.   Subsequently,  the  lift  and  drag 
coefficients  are  determined  through  the  use  of  Eq.  (15). 
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Fig.  17.   Flowchart  cf  computation  sequence. 
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The  normal  sequence  of  the  calculations  requires  a  test  for  the 
conditions  to  determine  whether  the  vortex  sheet  should  be  cut  and  the 
attached  vortex  spiral  be  shed.  Because  of  its  importance  and  the  need 
for  special  handling,  it  is  discussed  separately. 

In  discussing  the  problems  associated  with  the  use  of  discrete  vor- 
tices, it  has  been  pointed  out  that  the  elemental  vortices  must  be  placed 
at  the  center  of  the  sheet  segment  which  they  are  supposed  to  represent. 
This  led  to  the  method  of  red iscreti zation.  Thus,  following  the  test  for 
vortex  shedding  and  if  necessary,  after  the  cutting  of  the  vortex  sheet, 
all  attached  and  detached  sheets  are  red i scret i zed.   In  doing  so,  each  sheet 
is  separately  red iscreti zed  (see  Fig.  18).  The  discretization  of  the  de- 
tached sheets  is  handled  in  accordance  with  the  steps  described  in  Chapter 
II,  Section  E.   The  only  exception  to  the  said  process  is  in  the  treatment 
of  the  attached  sheets  near  the  nascent  vortex.   It  has  been  repeatedly 
pointed  out  that  the  conditions  near  the  separation  point  must  be  disturbed 
as  little  as  possible  by  redi screti zation  or  any  other  flow  of  events.   In 
accordance  with  this  principle,  the  nascent  vortex  and  its  immediate 
neighbor  are  not  allowed  to  change  their  strengths  or  positions. 

The  next  step  in  the  sequence  of  calculations  calls  for  circulation 
reduction  and  coalescence.   It  is  simply  hypothesized  that  every  vortex  in 
fhe  flow  field  loses  its  strength  in  an  amount  proportional  to  its  current 
strength.  This  hypothesis  is  based  on  heuristic  reasoning,  observations, 
and  numerous  calculations.   Further  justification  of  the  hypothesis  is  made 
later.  Here  only  the  actual  process  incorporated  into  fhe  computer  code  is 
descri  bed. 

During  each  computation  cycle  for  t  >  5,  the  strength  of  every  vor- 
tex, except  the  nascent  vortices,  is  reduced  according  to  the  following 
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expression 


T  (new  value)  =  (l-p/100)  V    (old  value) 
n  (       n 


(39) 


where  p  depends  on  the  location  x  =  Re{z  }  as  depicted  in  Fig.  19.  For 
r   r  n       n       r  3 

t  <  5,  p  is  varied  linearly  from  zero  to  its  final  value  of  unity. 
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9.   Variation  of  p  with  distance  for  t  >  5, 


As  discussed  in  connection  with  various  circulation  reduction  mech- 
anisms, there  are  two  other  processes  whereby  the  circulation  is  reduced. 
The  first  is  the  reduction  resulting  from  the  proximity  of  vortices  to  the 
cylinder.   Whenever  a  point  vortex  comes  closer  to  the  cylinder  than  a 
radial  distance  smaller  than  Ar  =  0.04,  that  point  vortex  is  assumed  to  be 
ann i h i I ated  and  removed  from  the  flow  field.   The  part  of  the  sheet  con- 
nected to  the  nascent  vortex  is  excluded  from  this  process.   The  annihila- 
tion of  circulation  in  this  manner  results  only  in  a  small  reduction  in  the 
total  circulation  and  cannot  account  for  a  40  to  50  percent  reduction  in 
the  strength  of  the  vortices  in  the  stable  region.   The  previous  applica- 
tions of  the  DVM  which  considered  this  annihilation  mechanism  together  with 
that  resulting  from  the  occasional  proximity  of  oppositely-signed  vortices 
(rediscretization  not  used)  obtained  only  a  mere  15  percent  reduction  in 
circulation  C67,  68],   This  fact  was  the  first  indication  of  the  need  for 
an  additional  circulation  reduction  mechanism. 
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The  second  of  the  two  mechanisms  concerns  the  entrainment  of  the 
tail  of  the  detached  sheet  into  the  vortices  across  the  wake.   For  the  sake 
of  discussion,  assume  that  an  attached  vortex  sheet  is  severed  at  a  suit- 
able point.  The  details  of  how  the  cutting  is  accomplished  is  discussed 
later.  The  tail  of  the  detached  vortex  sheet  is  drawn  into  the  region  be- 
tween two  attached  sheets  because  of  induction.   If  left  alone  that  tail 
eventually  develops  complex  kinks  and  spirals  and  makes  the  rediscreti za- 
tion  process  unnecessarily  complex.  As  proposed  by  Gerrard  C63],  the  cir- 
culation in  the  tail  gradually  ends  up  in  the  attached  sheet  across  the 
wake.   In  the  model  a  similar  procedure  is  employed.   That  portion  of  the 
tail  which  penetrates  into  the  mixing  region  beyond  a  line  tangent  to  the 
two  attached  vortex  spirals  (upstream  of  line  AA'  shown  in  Fig.  20)  is  re- 
moved from  the  fiow  field,  and  its  circulation  is  uniformly  subtracted 
from  that  of  the  sheet  across  the  wake.  Calculations  which  keep  track  of 
this  cross-wake  circulation  reduction  show  that  approximately  five  percent 
of  the  circulation  of  the  detached  sheet  is  transferred  across  the  wake 
by  this  mechanism.   It  should  also  be  noted  that  the  position  of  the  line 
severing  the  tail  of  the  detached  sheet  is  not  of  particular  significance. 
One  may  have  used  other  suitable  and  defensible  procedures  to  accomplish 
the  same  purpose. 

2.   Initiation  of  Asymmetry  and  the  Vortex  Shedding  Processes 

The  first  of  these  two  events  occurs  only  once  and  the  second  every 
time  period  a  vortex  is  shed.   As  such,  they  constitute  special  events  in 
the  computational  sequence  and  are   discussed  separately. 

a.   Asymmetry  Introduction 

In  nature  the  inception  of  asymmetry  and  the  shedding  of  vor- 
tices are  intrinsic  characteristics  of  flow  about  bluff  bodies.   The  causes 
of  asymmetry  and  the  subsequent  vortex  shedding  are  closely  linked  to  the 
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instabilities  in  the  shear  layers.  These  instabilities  may  be  brought  about 
by  the  slight  nonuni formities  in  the  ambient  flow,  randomly  distributed  tur- 
bulence, asymmetry  of  the  body,  etc.   It  appears  that  the  symmetric  state 
of  separated  flow  is  not  its  most  stable  state. 

Whether  the  flow  analysis  is  by  the  DVM,  finite  difference, 
finite  element  or  whatever,  asymmetry  must  be  artificially  introduced. 
Otherwise,  the  flow  field  remains  forever  perfectly  symmetric  and  the  vortex 
shedding  process  never  begins.   Figure  21  illustrates  the  results  of  the 
present  model  when  asymmetry  is  never  forced  to  develop. 

Although  the  onset  of  asymmetry  and  the  subsequent  evolution 
of  the  wake  have  been  observed  and  recorded  by  many  researchers,  little  is 
understood  about  the  controlling  mechanism.  The  question  of  whether  the 
asymmetry  begins  in  the  wake  and  is  fed  back  to  the  boundary  layer,  or 
vice  versa,  is  unknown.   Perhaps  both  occur  simultaneously.   Regardless, 
the  process  is  not  well  enough  understood  to  suggest  any  preferable  way 
to  artificially  initiate  its  occurrence.  Consequently,  the  method  utilized 
in  the  present  model  is  a  result  of  heuristic  reasoning  based  on  well 
known  experimental  observations  and  numerous  numerical  experiments  with 
the  model.   Before  describing  the  details  and  results  of  the  various  at- 
tempts, it  is  worthwhile  at  this  point  to  relate  some  of  the  general  ob- 
servations, or  lessons  learned,  from  the  investigation.   Firstly,  whatever 
method  of  asymmetry  introduction  is  chosen,  the  steady  state  results  of 
the  model  are  identical;  i.e.,  the  model  has  no  memory  of  the  conditions 
which  initiated  the  shedding  of  the  first  vortex.   Secondly,  as  stated 
earlier,  asymmetry  need  be  introduced  only  for  the  shedding  of  the  first 
vortex.   Following  the  shedding  of  the  first  vortex,  the  shedding  process 
continues  under  its  own  impetus.   Thirdly,  every  asymmetry  mechanism  in- 
volves a  trade-off  between  the  intensity  and  the  duration  of  the  perturbation 
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Fig.  21.   Evolution  of  flow  without  asymmetry  introduction. 
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applied.  Associated  with  the  question  of  duration  is  the  decision  of  when 
to  initiate  the  mechanism.   Fourthly,  there  is  an  optimum  time  during 
which  the  symmetric  flow  is  most  susceptible  to  the  onset  of  asymmetry. 
If  asymmetry  is  introduced  during  this  period,  the  intensity  and  duration 
of  its  appl ication  can  be  minimized.   Final ly,  of  the  unl  imited  number  of 
mechanisms  that  could  be  devised,  it  is  best  to  choose  the  simplest  method 
that  produces  the  desired  asymmetry. 

The  interactive  nature  of  the  vorticity  generation  mechanism 
and  the  vorticity  distribution  in  the  wake  suggest  that  asymmetry  might  be 
introduced  either  into  the  generation  mechanism  or  into  the  wake  itself. 
Methods  applied  to  the  vorticity  generation  method  involve,  for  example, 
artificially  increasing  the  computed  vorticity  flux  at  one  separation  point 
(say  by  10  percent)  and  simultaneously  decreasing  the  flux  at  the  other 
separation  point.   However,  sample  calculations  with  this  procedure  show 
that  the  method  is  not  effective.   Furthermore,  the  large  asymmetry  causes 
the  computed  separation  angles  to  undergo  unrealistic  excursions  from  their 
mean  positions. 

Introduction  of  asymmetry  into  the  wake  can  be  accomplished  by 
either  artificially  altering  the  strength  and/or  position  of  the  point  vor- 
tices on  the  existing  sheets.   Changing  the  strengths  of  the  vortices  in 
one  sheet  relative  to  the  other  results  in  a  sluggish  response  of  the  wake 
similar  to  that  described  above.  The  general  ineffectiveness  of  altering 
vortex  strengths  to  produce  the  downstream  convection  of  one  sheet  relative 
to  the  other  suggests  the  more  direct  approach  used  in  the  present  model. 
The  rationale  Is  to  initiate  the  downstream  convection  of  one  sheet  by 
actually  displacing  it  a  small  amount  each  computation  cycle  over  a  short 
period  of  time  so  that  when  the  mechanism  is  discontinued  the  sheet  con- 
tinues to  move  without  external  influence. 
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Numerous  numerical  experiments  have  shown  that  a  gradual  appli- 
cation and  buildup  of  the  displacement  increment  serves  to  provide  a  smooth 
transition  from  symmetric  to  asymmetric  conditions.  The  initiation  time 
and  duration  of  application  are  based  primarily  on  two  considerations. 
Firstly,  this  time  interval  should  correspond  to  the  period  during  which 
the  flow  is  most  susceptible  to  asymmetry.  Secondly,  the  initiation  of  the 
mechanism  should  begin  after  the  initial  peak  in  the  drag  curve.   Since  the 
cause  of  this  drag  overshoot  is  the  symmetric  growth  of  the  vortices,  the 
asymmetry  mechanism  should  not  be  allowed  to  influence  this  symmetric  de- 
velopment. 

In  the  following  the  specifics  of  the  asymmetry  introduction 
method  are  presented.   During  each  computation  cycle  within  the  interval 
5.0  <_  t  <_  9.0  the  vortices  of  the  upper  attached  sheet  are  displaced  down- 
stream an  amount  given  by  (see  Fig.  22) 


Ax  =  0.08  {1.0  -  cos[ir(t-5)/4]} 


(40) 


0.16  - 


0.08  - 


0.00 


Fig.  22.   Variation  of  the  asymmetry  perturbation  Ax  with  time. 
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Several  points  are  worthy  of  discussion.   Firstly,  the  displacement  given 
by  Eq.  (40)  is  not  unique  and  any  other  form  of  displacement  function  could 
have  accomplished  the  same  purpose.  The  particular  choice,  depicted  by 
Fig.  22,  introduces  the  disturbances  gradually  and  continuously  without 
shocking  the  flow  field.   Lastly,  the  total  displacement  applied  to  the 
sheet  during  the  time  interval  5  <_  t  <_  9  is  only  32  percent  of  the  radius 
of  the  cy I inder. 

b.   Vortex  Shedding  Process 

Vortex  shedding  is  understood  as  the  mechanism  whereby  the  feed- 
ing layer  is  cut.   The  specific  details  of  the  conditions  which  dictate  the 
separation  of  a  spiralling  sheet  from  its  feeding  layer  are  not  known  and 
have  yet  to  be  determined  experimentally.  The  numerical  model  is  not  in 
need  of  such  a  cut  were  it  not  due  to  the  limitations  of  computation  time. 
If  the  flow  field  is  allowed  to  evolve  asymmetrically,  the  flow  pattern 
shown  in  Fig.  23  results.   This  figure  shows  not  only  the  evolution  of  the 
doubly-connected  vortex  sheets  but  also  the  drawing-in  of  the  tail  of  the 
shedding  vortex  into  the  mixing  region.  The  continuation  of  the  calculations 
in  this  manner  are  prohibitive  even  with  the  largest  computers.   However, 
it  is  instructive  in  two  ways  to  examine  the  evolution  of  the  wake  in  this 
manner.   Firstly,  the  uncut  sheet  begins  to  change  the  sign  of  its  curva- 
ture as  the  lift  force  goes  through  its  extremal  (see  point  A  of  Fig.  23). 
The  change  of  curvature  is  in  fact  the  beginning  of  the  drawing-in  process 
of  the  shear  layer  across  the  wake.   Subsequently,  a  strong  inflection  and 
kink  develop  drawing  the  sheet  closer  into  the  formation  region.   Simul- 
taneously, the  lift  force  increases  towards  its  next  extremal.   The  gradual 
nature  of  this  process  does  not  permit  one  to  identify  a  specific  moment 
at  which  the  shedding  should  commence.   Thus,  the  examination  of  other  flow 
features  exhibiting  stronger  reaction  to  the  precise  moment  of  shedding  are 

88 


Fig.  23.   Evolution  of  flow  without  cutting  the  vortex  sheet. 
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necessary.  Numerous  calculations  with  symmetric  and  asymmetric  evolution 
of  the  wake  and  with  and  without  the  rediscreti zation  of  the  shear  layers 
have  shown  that  the  sheet  to  be  cut  bears  its  minimum  circulation  per  unit 
length  at  a  time  immediately  after  dT/dt  goes  through  a  minimum  (see  Fig. 
24).   In  other  words,  the  connecting  sheet  is  weakest  and,  therefore, 
most  susceptible  to  breaking  when  its  circulation  is  minimum.   It  is  on 
the  basis  of  this  fact  that  the  sheet  is  cut  immediately  after  the  rate 
of  circulation  reaches  its  minimum.   At  the  time  of  cut,  the  lift  force 
is  about  50  percent  of  its  preceding  extremal  (see  Fig.  25). 

The  actual  cutting  process  is  accomplished  by  removing  a  single 
point  vortex,  for  only  one  time  step,  from  the  sheet  at  a  distance  s  =  0.4 
(Fig.  26).   Subsequent  calculations  separately  rediscretize  the  detached 
sheet  and  the  part  sti I  I  connected  to  the  nascent  vortex.   As  the  flow 
evolves,  the  attached  sheet  continues  to  grow  rapidly  to  form  a  new  spiral 
while  the  detached  sheet  moves  further  downstream.   At  each  time  step  the 
computer  code  tests  the  condition  as  to  whether  the  new  connected  sheet 
should  be  cut.   At  the  time  the  said  condition  is  satisfied,  the  vortex 
sheet  is  cut  as  described  above.   It  is  clear  that  the  shedding  of  the  vor- 
tices and  the  cutting  process  alternate  with  a  time  period  equal  to  half 
the  Strouha I  period.  When  a  vortex  sheet  is  cut,  say  at  the  lower  side  of 
the  cylinder,  the  sheet  which  was  previously  detached  from  the  top  part  of 
the  cylinder  is  coa I esced  into  a  single  vortex  at  its  center  of  gravity. 
At  the  time  of  its  coalescence,  the  core  of  the  vortex  spiral  has  moved  a 
distance  at  least  seven  radii  downstream.   The  calculations  have  conclu- 
sively shown  that  the  process  of  coalescence  at  this  distance  from  the 
cy I  inder  does  not  a  Iter  the  continu  ity  of  any  variable  (e.g.,  lift,  drag, 
circulation,  etc.). 
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Fia.  25.   Variation  of  lift  during  a  cutting  cycle. 
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Fig.  26.   Cutting  of  a  vortex  sheet. 


93 


B.   DISCUSSION  OF  RESULTS  AND  SENSITIVITY  ANALYSIS 

In  the  following,  the  results  obtained  with  the  numerical  model  (re- 
ferred to  hereafter  as  the  standard  run)  are  presented.   Following  this, 
the  variation  of  some  of  the  important  kinematic  and  dynamic  parameters 
to  the  variation  in  others  is  discussed.  The  purpose  of  these  discussions 
is  to  partly  establish  the  stability  of  the  results,  partly  to  relax  the 
range  of  variability  of  the  nond isposab le  parameters,  and  partly  to  estab- 
lish the  cause  and  effect  relationships.  At  this  time,  no  distinction 
is  offered  concerning  the  choice  of  parameters  in  the  model  and  the  calcu- 
lated results.   Furthermore,  no  distinction  is  made  between  prediction  and 
inspired  hindsight. 

I .   Results  Obtained  with  the  Standard  Run 

The  standard  run  is  computed  using  At  =  0.125,  p  =  1.0  (see  Fig.  19) 
and  with  the  asymmetry  introduction  and  the  shedding  mechanisms  previously 
descr i  bed. 

The  evolution  of  the  wake  over  one  arbitrarily  chosen  cycle  is  shown 
in  Figs.  27  and  28.   Figure  27  depicts  the  near  wake  region  and  Fig.  28 
shows  the  far  wake  vortex  street. 

The  variations  of  the  lift  and  drag  coefficients  with  t  are  shown 
in  Figs.  29  and  30.  The  drag  coefficient  reaches  a  maximum  at  t  -   4.2  and 
then  rapidly  decreases  to  its  ultimate  mean  value  of  about  1.2.   The  fluc- 
tuations of  CR  have  a  frequency  twice  the  vortex  shedding  frequency  and  an 
amplitude  of  about  0.06.   The  overshoot  of  the  drag  coefficient  is  a  direct 
consequence  of  the  rapid  accumulation  of  vorticity  in  the  symmetrically 
growing  vortex  spirals.   Following  the  shedding  of  the  first  vortex,  Cn 
drops  sharply  and  the  lift  coefficient  increases,  ultimately  reaching  a  mean 
amplitude  of  about  0.65.   The  slight  variations  in  the  oscillations  of  Cn 
and  C  come  from  various  sources.   The  vortices  absorbed  during  the  cross- 
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Fig.  27.   Evolution  of  flow  in  the  near  wake. 
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Fiq.  27  (con't).   Evolution  of  flow  in  the  near  wake. 
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Fig.  27  (con't).   Evolution  of  flow  in  the  near  wake. 
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wake  cancellation  mechanism  may  be  one  more  or  one  less  during  two  different 
cycles.  The  increase  of  theotherwise  finite  length  of  the  wake  does  not 
permit  the  vortex  sheets  to  follow  exactly  identical  paths  from  one  cycle 
to  another.  Regardless,  the  said  perturbations  in  Cn  and  C.  are  quite 
small  and  well  within  the  accuracy  of  the  numerical  model. 

The  amplitude  of  the  lift  coefficient  is  certainly  within  the 
range  of  lift  coefficients  reported  in  the  literature.  However,  it  is  not 
possible  to  find  a  truly  two-dimensional  laboratory  experiment  to  compare 
the  calculated  results  with  greater  def i n i ti veness.   More  is  said  later 
about  the  dependence  of  the  amplitude  of  C.  on  other  parameters,  in  par- 
ticular on  the  circulation  reduction  mechanism. 

I       Figure  31  shows  a  comparison  of  the  calculated  and  measured  drag 
coefficients.   The  results  are  strikingly  similar  and  may  be  taken  as  a 
first  indication  of  the  workings  of  the  model. 
The  oscillations  of  the  separation  points  in  terms  of  the  angles 
»  measured  from  9=0  (not  necessarily  the  instantaneous  position  of  the  stag- 
nation point)  are  shown  in  Fig.  32.   The  separation  angle  for  a  given  sepa- 
ration po  int  osci I  I ates  near  I y  s  inuso  i  da  I  I y  with  an  amp  I  itude  of  A6  =  3 
degrees  about  a  mean  va I ue  of  9  -   11   degrees.   The  frequency  of  oscilla- 
tion of  9   is  identical  to  that  of  the  shedding  of  vortices, 
s  3 

The  stagnation  point  does  not  remain  at  9  =  0  even  in  steady  flow 
about  a  stationary  cylinder.   Figure  33  shows  the  oscillation  of  the  stagna- 
tion point.   Similar  oscillations  have  been  reported  by  Dwyer  and  McCroskey 
\J\~\.      The  separation  angles  relative  to  the  instantaneous  position  of  the 
stagnation  point  are  shown  in  Fig.  34.   A  comparison  of  Figs.  32  and  33 
show  that  the  oscillations  of  fhe  separation  and  stagnation  points  are  180 
degrees  out  of  phase.   In  other  words, J f  the  stagnation  pornt  moves  up, 
say  to  9  =  2,  the  separation  point  at  the  upper  half  cf  the  cylinder  moves 
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backwards  from  a  mean  position  9  =  77  to  8  =74  degrees.  This  fact  is 
also  in  conformity  with  previously  reported  results  C7lH. 

The  rate  of  change  of  circulation  at  both  the  upper  and  lower  sepa- 
ration points  is  shown  in  Fig.  35.  Because  the  generation  and  shedding  of 
circulation  are  the  essence  of  the  entire  bluff  body  problem,  dT/dt  has  a 
special  significance  among  all  other  results  predicted  by  the  numerical 
model.  A  cursory  examination  of  this  figure  shows  that  there  is  a  perfect 
similarity  of  the  rates  of  circulation  of  the  upper  and  lower  nascent  vor- 
tices.  Secondly,  the  variations  in  dF/dt  from  one  cycle  to  another  as  t 
increases  from  40  to  200  are  very  minor,  indicating  that  the  vortex  street 
may  be  regarded  as  infinitely  long  for  all  intents  and  purposes  for  t 
larger  than  about  40. 

The  mean  value  of  |dT/dt|  is  about  1.2.   Noting  that  the  mean  base 

pressure  C  ,  is  related  to  IdT/dtl  by  [72] 

pb  '     '   ' 

C  ,  =  1.0-2  Idf/dtl  (41  ) 

pb  '     ' 

The  present  calculations  yielded  C  .  =  -1.35  which  compares  favorably  with 
that  reported  by  Roshko  DG  at  high  subcritical  Reynolds  numbers  [C   , 
(experimental)  =  -1.36  for  20,000  <  Re  <  100,000]. 

The  amplitude  of  oscillation  of  dl/dt  is  about  0.18  which  indicates 
that  the  velocity  at  the  separation  point  varies  from  1.41  to  1.65. 

All  of  the  parameters  cited  above, namely  lift  and  drag  coefficients, 
separation  and  stagnation  angles,  and  dlYdt  oscillate  with  a  frequency 
identical  to  the  vortex  shedding  frequency  f.   The  Strouhal  number,  St  = 
fyD/U  is  found  to  be  St  -   0.205.   The  experimental  values  of  the  Strouhal 
number  are  reported  to  lie  between  0.195  and  0.210  for  the  Reynolds  number 
range  from  2,000  to  400,000  (see  Fig.  36).   Thus,  the  predictions  of  the 
model  are  in  conformity  with  the  measurements  in  the  range  of  Reynolds 
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Fig.  36.   Strouhal  number  versus  Reynolds  number  H73]. 
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numbers  where  the  DVM  is  applicable.  As  is  noted  later,  the  Strouhal  number 
is  the  least  sensitive  of  all  the  parameters  to  the  variations  of  the  non- 
disposable  parameters. 

One  of  the  often  discussed  characteristics  of  the  vortex  street, 
particularly  in  the  stable  region,  is  the  relative  spacing  of  the  vortices. 
Von  Karman's  classic  stability  analysis  gave  h  /£  =  0.281  (see  Fig.  37). 
Since  von  Karman,  the  said  ratio  has  been  determined  experimentally  by  a 

great  number  of  investigators  and  it  has  been  found  that  h  /£  varies  from 
3  a  v  v 

about  0. 19  to  0.3. 

The  lateral  and  transverse  vortex  spacings  obtained  with  the 
standard  run  are  shown  in  Fig.  37.   The  asymptotic  values  of  I     and  h  are 
well  within  the  experimentally  determined  values  C74H.   In  the  stable  region 
of  the  street,  the  lateral  spacing  decreases  somewhat  and  the  transverse 
spacing  gradually  increases  from  about  0.63  to  its  asymptotic  value  of  about 
1.05.   The  ratio  h  It      in  the  asymptotic  range  is  seen  to  be  0.23,  a  value 
which  compares  well  with  those  reported  experimentally  f.74].   Evidently, 
it  is  not  meaningful  to  think  of  vortex  spacings  within  the  formation 
reg  ion . 

All  discrete  vortex  models  of  separated  flow  about  bluff  bodies 
indicate  the  need  for  greater  loss  of  circulation  in  the  wake  in  order  to 
bring  a  closer  correspondence  between  the  calculated  and  measured  separa- 
tion points,  lift  and  drag  coefficients,  dr/dt,  etc.   For  example,  Deffen- 
baugh  and  Marshall  C47]  obtained  a  mean  separation  angle  of  67  degrees  and 
a  mean  drag  coef f ic ient  of  0.9  even  though  the  smaller  separation  angle 
should  have  yielded  larger  drag  than  that  found  experimentally  due  to  the 
increased  wake  size.   In  the  present  model  all  conceivable  circulation 
reduction  mechanisms  have  been  taken  into  consideration.   Their  variations 
with  x  or  t  may  not  be  unique  or  in  conformity  with  those  obtained 
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experimentally  but  their  characteristics  follow  all  vortex  shedding  mech- 
anisms so  far  advanced.   Figure  38  shows  the  circulation  retained  by  a 
vortex  sheet  from  its  inception  to  its  subsequent  coalescence  and  convec- 
tion into  the  far  wake.  Also  shown  in  this  figure  are  the  two  experimen- 
tal points  reported  by  Bloor  and  Gerrard  [69]  and  the  mean  curve  drawn 
through  the  data  reported  by  Schmidt  and  Tilmann  [JO'2   (see  Fig.  15).  The 

I  comparison  of  a  I  I  the  results  is  better  than  expected  in  view  of  the  fact 
that  the  experimental  determination  of  circulation  through  direct  or 
indirect  methods  is  an  extremely  difficult  and  approximate  task.   In  par- 
ticular, the  ultrasonic  sound  beam  technique  used  by  Schmidt  and  Tilmann 
yields  the  vortex  strengths  only  approximately  in  the  stable  and  unstable 
region  of  the  vortex  street.   There  are  no  measurements  of  the  vortex 
strengths  in  the  formation  region.   The  foregoing  does  not  imply  that  the 
particular  variation  assigned  to  the  circulation  loss  mechanism  (see  Fig. 
19)  is  unique.   In  fact,  the  only  major  nond isposab le  parameter  in  the 
entire  numerical  model  is  the  form  assigned  to  the  circulation  loss.   The 
need  for  circulation  loss  has  been  amply  demonstrated.   The  form  of  the 
particular  relationship  can  only  be  adjusted  with  hindsight  by  comparing 
all  of  the  predictions  of  the  model  with  those  obtained  experimentally. 
Because  of  the  fact  that  the  differences  in  measured  and  calculated  circu- 
lations are  further  downstream  in  the  wake,  the  use  of  the  experimentally 
suggested  p  variation  does  not  in  any  way  affect  lift,  drag,  separation 
angles,  dT/dt,  and  the  Strouhal  number.   This  wil I  become  more  evident  in 
the  course  of  the  discussion  of  the  sensitivity  of  the  model.   Calculations 
nof  reported  herein  with  more  exotic  forms  of  p  have  ascertained  the 
validity  of  the  insignificance  of  variation  of  p  beyond  the  stable  region. 
2.   Sensitivity  Analysis 

This  section  deals  with  the  effect  of  the  variation  of  one  parameter 
abouts  its  value  used  in  the  standard  run  on  the  remaining  major  parameters. 
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All  numerical  analyses  must  use  a  finite  mesh  size  or  time  interval. 
The  magnitude  of  these  quantities  might  affect  the  results  in  various  ways. 
In  the  DVM  a  large  At  might  bring  in  curvature  effects  particularly  in  flow 
regions  where  the  radius  of  curvature  of  the  spiral  is  small.  This  effect 
is  brought  about  by  the  convection  of  the  vortices  along  straight  lines  co- 
incident with  their  instantaneous  velocities.  Thus  At  should  be  kept  reason- 
ably small  to  allow  the  vortices  to  follow  their  natural  paths.  A  very 
small  At,  on  the  other  hand,  requires  prohibitively  large  computation  times 
and  defeats  the  purpose  of  the  numerical  experiments. 

The  standard  run  was  made  with  At  =  0.125.  Additional  runs  were 
made  varying  At  from  0.1  to  0.2.   The  resulting  drag  and  lift  coefficients 
are  presented  in  Figs.  39  and  40.   The  differences  in  the  two  coefficients 
for  all  va I ues  are  neg I igi  b ly  sma II  and  certain  I y  within  the  range  of 
accuracy  expected  from  numerical  calculations.  Other  parameters  such  as 
separation  and  stagnation  angles,  circulation  loss,  Strouhal  number,  etc. 
remain  essentially  identical  to  their  corresponding  values  obtained  with 
At  =  0. 125. 

The  role  played  by  the  loss  of  circulation  in  both  laboratory  and 
numerical  experiments  has  already  been  discussed.  Here  attention  is  devoted 
to  the  selection  of  the  particular  value  of  p  and  the  effect  of  its  varia- 
tion on  the  remainder  of  the  predictions.   Several  hundred  hours  of  calcu- 
lations have  shown  that  one  can  arrive  at  an  approximate  value  of  p  by 
comparing  the  predicted  and  measured  values  of  the  vortex  strengths  and  the 
i f t  and  drag  coefficients.   None  of  these  parameters,  however,  exhibit  the 
same  degree  of  variation  with  p  at  corresponding  times  and  downstream  dis- 
tances.  For  example,  the  comparison  of  the  vortex  strengths  in  the  stable 
region,  where  most  of  the  measurements  are  made,  is  not  very  meaningful 
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Fig.  39.   Drag  coefficient  versus  time  for  At=0.l,  0.125,  0.15  and  0.2. 
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Fig.  40.   Lift  coefficient  versus  time  for  At=0.l,  0.125,  0.15  and  0.2. 
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since  the  effect  of  such  vortices  on  the  lift  and  drag  coefficients  is  very 
smal I.  Thus,  it  was  decided  to  choose  the  p  value  such  that  the  terminal 
value  of  the  drag  coefficient  conforms  to  that  measured  in  the  subcritical 
region.   No  other  experimental  data  was  used  in  the  analysis.   Numerical 
experiments  with  various  values  of  p  have  shown  that  p  =  1.0  yields  (L  - 
1.2.  Having  determined  a  value  of  p,  it  is  necessary  to  examine  its  varia- 
tion on  all  of  the  predictions  of  the  model. 

The  drag  coefficients  obtained  with  p  =  0,  p  =  0.5,  p  =  I,  and 
p  =  2  are  shown  in  Fig.  41.   The  insert  in  Fig.  41  shows  the  variation  of 
C  with  p  in  the  asymptotic  range  of  CL.   Evidently  a  100  percent  change  in 
p  about  p  =  I  caused  only  about  10  percent  change  in  C~.   Consequently,  the 
tuning  of  the  model  through  the  use  of  CL  =  1.2  is  not  as  critical  as  it 
may  have  appeared  at  first  sight.   Figure  41  also  shows  that  a  circulation 
loss  mechanism  is  necessary  (see  curve  for  p  =  0)  if  the  numerical  model 
is  to  yield  results  in  conformity  with  the  experiments.   As  noted  earlier, 
the  need  for  circulation  loss  is  also  evidenced  by  the  vortex-strength 
measurements  (see  Fig.  38).   Thus,  the  important  question  in  the  evolution 
of  the  model  is  not  whether  there  should  be  a  circulation  loss  but  rather 
what  should  the  mechanism  and  magnitude  of  this  loss  be. 

The  variation  of  the  lift  coefficient  with  p  is  shown  in  Fig.  42. 
It  is  immediately  apparent  that  the  lift  coefficient  is  far  more  sensitive 
to  p  than  C_.   Regardless  of  its  sensitivity,  however,  C.  remains  well 
within  the  experimentally  reported  values.  One  could  fine-tune  the  value 
of  p  by  matching  the  measured  and  calculated  values  of  C  .   Unfortunately, 
thi.s  is  not  possible  since  there  are  no  reliable  C  values,  particularly  in 
the  region  of  high  subcritical  Reynolds  numbers.   Consequently,  no  attempt 
is  made  to  slightly  increase  or  decrease  p  from  p  =  I  to  obtain  a  specific 
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Fig.  41.   Drag  coefficient  versus  time  for  p=0.0,  0.5,  1.0  and  2.0. 
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Fig.  42.   Lift  coefficient  versus  time  for  p=0.0,  0.5,  1.0  and  2.0. 
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C.  value.  The  importance  of  Fig.  42,  rather,  lies  in  the  fact  that  th< 


numerical  model  clearly  shows  the  sensitivity  of  C  to  circulation  loss  and 
offers  an  indirect  explanation  of  the  difficulties  encountered  in  the 
laboratory  experiments.   The  findings  of  the  model  can  also  be  taken  as  a 
means  of  reducing  the  transverse  force  by  artificially  Increasing  the 
circulation  dissipation  in  vortices  (e.g.,  use  of  fairings,  wire  screens, 
etc . ) . 

The  Strouhal  number  displays  a  relatively  small  variation  wifh  p 
(see   Fig.  43).   This  explains  the  good  correspondence  obtained  between 
measured  and  calculated  values  of  St  even  when  the  circulation  reduction 
does  not  exceed  15  percent  (Clements  C67H,  Sarpkaya  C 6 8J ) % 


0.31-  St 


0.2 


0.1 


0.0  1.0 

Fig.  43.   Variation  of  the  Strouhal  number  with  p. 
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It  is  of  interest  to  point  out  the  relationship  between  the  strength 
and  motion  of  vortices  in  the  near  wake  and  the  small  variations  in  the 
Strouhal  number.   When  the  detached  vortex  is  relatively  stronger  (p  <  I), 
then  it  lingers  longer  near  the  cylinder  by  virtue  of  the  velocities  in- 
duced on  it  by  the  equally-strong  vortex  across  the  wake.   This  in  turn  in- 
creases the  period  of  vortex  shedding  or  decreases  the  Strouhal  number. 
For  weaker  vortices  (p  >  I)  the  inverse  is  true. 


In  the  following,  the  sensitivity  of  the  lift  and  drag  coefficients 
and  the  Strouhal  number  to  the  variations  in  the  amplitude  of  p  has  been 
discussed.   In  doing  so,  the  shape  of  the  dissipation  function  has  not 
been  changed.   From  a  scientific  point  of  view,  the  next  to  consider  is 
the  effect  of  the  variation  of  the  dissipation  function  on  the  major  predic- 
tions of  the  model.   For  this  purpose  an  extensive  study  was  undertaken  to 
acquire  a  deeper  understanding  of  the  dissipation  particularly  in  the 
near  wake  region  of  the  cylinder.   In  doing  so,  not  only  the  terminal 
value  of  the  drag  coefficient  but  also  other  parameters  such  as  the  lift 
coefficient,  separation  points,  and  the  amount  of  circulation  retained  by 
the  concentrated  vortices  were  carefully  observed.   Here  only  the  high- 
lights of  these  observations  can  be  described.   The  use  of  a  dissipation 
function  such  as  p  =  a/x   (for  x  >  I)  where  a  is  a  constant  shifts  the  re- 
gion of  larger  dissipation  closer  to  the  rear  of  the  cylinder.   For  example, 
the  results  obtained  with  a  =   2.0  have  shown  that  the  lift  and  drag  coeffi- 
cients show  no  more  than  a  10  percent  variation  from  those  obtained  with 
the  standard  model.   The  separation  points  remained  at  almost  exactly  the 
same  average  values.   The  percent  circulation  retained  by  the  concentrated 
vortices  increased  for  obvious  reasons.   In  fact,  the  relative  vortex 
strengths  at  comparable  distances  were  somewhat  larger  than  those  measured 
by  Gerrard  [24^.   The  Strouhal  number  decreased  from  0.205  to  0.193,  in- 
dicating that  not  only  the  vortices  in  the  formation  region  but  also  those 
further  downstream  determine  the  transport  velocity  of  the  vortices  and 
hence  the  Strouhal  number. 

The  experience  gained  from  this  exploratory  investigation  has  shown 
"hat  there  should  be  a  fairly  large  and  uniform  dissipation  in  the  near 
wake,  say  within  the  formation  region,  and  then  decrease  gradually  with 
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distance.  However,  without  careful  measurements  of  the  vortex  strengths 
as  a  function  of  time  or  distance  it  may  never  be  possible  to  determine 
the  precise  shape  of  the  dissipation  function. 

To  further  determine  the  sensitivity  of  the  model  to  dissipation  and 
the  relationship  between  some  of  the  most  important  parameters  (e.g., 
separation  angles,  vortex  strengths,  and  Strouha I  number),  an  additional 
series  of  numerical  experiments  have  been  performed.  Specifically,  the 
separation  angle  was  determined  not  in  the  manner  described  in  the  stan- 
dard model  but  rather  by  suitably  extrapolating  the  velocity  profiles 
beyond  the  point  of  maximum  velocity.   The  use  of  such  an  extrapolation 
alleviated  the  influence  of  the  nascent-vortex  proximity  effect  on  the 
separation  point  and  resulted  in  larger  separation  angles  (6  -   81  degrees). 
When  the  shape  of  the  dissipation  function  was  kept  exactly  as  in  the 
standard  run,  the  drag  coefficient  reached  an  ultimate  value  of  about  1.3 
and  the  lift  coefficient  a  value  of  about  0.8.   The  Strouhal  number  de- 
creased to  about  0.19.   It  was  thus  apparent  that  the  value  of  p  in  the 
standard  dissipation  should  be  somewhat  increased,  say  to  p  -  1.3.   Such 
an  increase  in  p  would  decrease  both  C  and  C.  and  increase  the  Strouhal 
number  without  affecting  the  separation  and  stagnation  points.   This 
numerical  experiment  suggests  that  it  is  possible  to  match  the  observed 
and  calculated  separation  points  and  resistance  to  fluid  motion  with  one 
nond i sposab le  parameter.   No  attempt  has  been  made  to  find  a  precise  value 
or  shape  for  the  dissipation  function  because  of  the  primary  concern  with 
the  interrelationships  between  various  parameters.   The  fact  has  been 
demonstrated  that  the  physical  mechanism  governing  the  flow  is  not  affected 
by  the  actual  value  of  p  or  Cp  and  C.  .   It  is  evident  from  the  foregoing 
that  the  numerical  modelling  of  a  turbulent  wake  does  not  require  a 
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dissipation  based  on  turbulent  diffusion  and  that  diffusion  could  be 
quantified  only  by  additional  measurements  which  will  complement  the 
numerical  experiments.  As  is  often  the  case  in  many  aspects  of  fluid 
mechanics,  the  question  seems  to  have  been  reduced  once  again  to  the 
understanding  of  the  physics  of  turbulence. 


V.   SUMMARY  OF  THE  RESULTS 

The  numerical  model  based  on  the  rediscreti zation  of  vortex  sheets  and 
the  mutual  interaction  between  the  boundary  layers  and  the  wake  yielded 
numerous  results  which  may  be  summarized  as  follows: 

1.  A  discrete  vortex  model  of  flow  past  a  stationary  cylinder  may  be  for- 
mulated with  only  one  nondisposab le  parameter.   Experiments  and  the  numeri- 
cal results  show  the  necessity  of  circulation  reduction  but  not  necessarily 
the  exact  mechanism  by  which  it  is  accomplished. 

2.  Based  on  a  postulated  circulation  reduction  mechanism,  all  predictions 
of  the  model  are  in  conformity  with  those  obtained  experimentally.  The 
fact  that  al I  and  not  only  a  few  of  the  parameters  agree  with  the  experi- 
mental measurements  adds  further  credibility  to  the  assumptions  made  and 
to  the  power  of  prediction  of  the  model. 

3.  A  concerted  effort  is  made  to  establish  a  balance  between  simplicity 
and  accuracy.   To  this  end,  simple  linear  relationships  have  been  preferred 
in  lieu  of  more  complex  ones. 

4.  One  of  the  most  important  aspects  of  the  numerical  model  is  its  flexi- 
bility.   Practically  all  parameters  have  a  wide  stable  range  and  allow 

the  user  to  perform  numerical  experiments  to  seek  a  physical  understanding 
of  the  intricate  relationships  between  two  or  more  parameters. 

5.  The  demonstration  of  the  applicability  of  the  discrete  vortex  model 
to  the  analysis  of  flow  about  a  bluff  body  with  mobile  separation  points 
shows  that  similar  models  can  be  constructed  with  ease  for  bodies  with 
fixed  separation  points.   This  allows  one  to  study  such  problems  as  vortex 
shedding  from  ship  bilge  keels,  stall  over  aircraft,  fuel  slushing  in  large 
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containers,  etc.,  in  naval  hydrodynamics  and  aeronautics. 
6.  The  model,  as  it  stands,  can  be  used  in  the  analysis  of  a  great  many 
flow-structure  interaction  problems.  Of  particular  interest  is  the  analysis 
of  a  circular  cylinder  undergoing  forced  or  self-excited  transverse  oscil- 
lations. Clearly,  the  numerical  experiments  with  oscillating  cylinders 
can  not  only  shed  additional  light  on  the  mechanisms  which  bring  about  such 
oscillations  but  also  increase  the  universality  of  the  model. 
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VI.  APPLICATION  OF  THE  DISCRETE  VORTEX  MODEL 
TO  SELF-EXCITED  TRANSVERSE  OSCILLATIONS 

A.   INTRODUCTION 

Numerous  experiments  have  shown  that  (see  Parkinson  [_"75j  for  an  ex- 
tensive review)  when  the  natural  frequency  of  a  bluff  body  is  close  to 
the  vortex-shedding  frequency,  the  body  may  undergo  self-excited  trans- 
verse oscillations.   The  primary  consequences  of  these  oscillations  are 
as  follows:   (i)   The  vortex  shedding  frequency  locks  on  to  the  natural 
frequency  of  the  body.   (ii)   The  amplitudes  of  oscillation  and  the 
transverse  force  increase  rapidly  up  to  the  point  of  perfect  synchroni- 
zation where  the  maximum  amplitude  is  attained.  Subsequently,  the  oscil- 
lations decrease  either  rapidly  or  abruptly.  The  abrupt  drop  in  amplitude 
is  a  consequence  of  hysteresis.  The  reasons  leading  to  the  hysteresis 
are  not  yet  known,   (iii)   The  transverse  force  leads  the  displacement 
by  a  phase  angle  which  increases  gradual ly  at  first  to  about  50  degrees 
and  then  rapidly  to  about  130  degrees,   (iv)   At  the  end  of  the  synchroni- 
zation region,  the  vortex  shedding  frequency  jumps  to  that  governed  by 
the  Strouhal  relationship. 

The  studies  concerned  with  the  understanding  and  prediction  of  fhe 
synchronization  phenomenon  may  be  classified  in  three  categories.   The 
first  and  most  extensive  category  has  been  the  experimental  observations 
and  measurements  with  forced  or  self-excited  oscillations  of  circular 
cylinders,  cables,  and  a  small  group  of  other  bluff  bodies.   The  second 
category  of  investigations  dea I f  primarily  with  the  analysis  of  the 
structural  aspects  of  the  equations  governing  the  vibrations  by  making 
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ad  hoc  assumptions  about  driving  fluid  forces.   Finally,  the  third  cate- 
gory concentrated  on  the  understanding  of  the  flow  field  and  the  nature 
of  the  fluid  forces.   It  has  not  yet  been  possible  to  understand  the 
coupled  structural  and  fluid  mechanical  mechanisms  to  explain  the  reasons 
leading  to  the  self-excited  oscillations. 

Mathematical  models  based  on  the  assumption  of  a  nonlinear  oscillator 
have  been  developed.  Among  those,  the  most  noteworthy  is  that  proposed 


by  Hartlen  and  Currie  [76H.  Their  model  predicts  some  of  the  observed 
features  of  the  vortex-induced  oscillations.   However,  the  model  requires 
a  number  of  nondisposab le  parameters  and  gives  no  clues  about  the  cause 
and  effect  relationships.   It  has  been  realized,  as  suggested  by 
Parkinson  C75H  among  others,  that  if  any  understanding  of  the  phenomenon 
is  to  be  achieved,  theoretical  developments  and  numerical  experiments 
must  be  compared  with  laboratory  experiments  to  guide  and  complement  each 
other.   It  is  with  this  suggestion  in  mind  that  the  discrete  vortex  model 
has  been  applied  to  the  prediction  of  the  characteristics  of  self-excited 
osci  I  I  at  ions  of  a  c  i  rcu I ar  cy I  i  nder.   I n  the  f ol  low  i  ng,  f i  rst  the 
description  of  the  mathematical  foundations  of  the  problem,  then  the 
specific  details  of  the  numerical  procedure,  and  finally  the  discussion 
of  the  results  and  a  summary  are  presented. 

B.   FORMULATION  OF  THE  PROBLEM 

The  equation  of  motion  for  an  e lastica I  I y-mounted  and  linearly-damped 
cylinder  of  unit  length  (see  Fig.  44)  may  be  written  in  dimensional  form 
as 

2        2 
my  +  2m  co  Cy  +  m  <jj   y  =  px  U  c  C,  (42) 

n        n  '    t      L 
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where  m  represents  the  mass  of  the  cylinder;  y,  the  displacement;  oj  ,  the 
natural  circular  frequency;  £,  the  damping  ratio;  p,,  the  fluid  density; 
c,  the  radius  of  the  cylinder;  and  C,  ,  the  instantaneous  value  of  the 
lift  coefficient.   In  Eq.  (42),  y  and  y  represent  respectively  the  first 
and  the  second  derivatives  of  y  with  respect  to  time. 
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Fig.  44.   E lastica I ly-mounted,  linearly-damped  cylinder. 


Defining  f   as  the  vortex  shedding  frequency  for  the  stationary 
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Eq.    (42)    may   be   written    in    normalized    form  as 


n  +  n  +   ( )      0  =   277  5     a  c, 

GO  go  o      o    L 

o  o 


(44) 
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Because  of  the  need  to  distinguish  the  stationary  and  oscillating  vortex 
shedding  frequencies,  and  hence  the  Strouhal  numbers,  the  previously  and 
commonly  used  notation  of  St  for  the  Strouhal  number  is  modified  as 
denoted  above. 

The  evaluation  of  Eq.  (44)  requires  the  instantaneous  value  of  the 
lift  coefficient.   In  the  present  model  it  is  evaluated  through  the  use 
of  the  complex  velocity  potential  modified  to  take  into  account  the  motion 
of  the  cylinder.  Thus,  one  has 
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w(z)  =  -[(z-z  )+  ->  ^-Vr  {£n(z-z  )-£n[(z-z  )- 

o   z-z    2tt  L^j      n       n        o   
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in  which  z  denotes  the  instantaneous  position  of  the  center  of  the  cylin- 
der and  z  ,  its  velocity  (see  Fig.  45).  The  use  of  Eq.  (45)  together  with 
the  generalized  Blasius  theorem  yields 
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Fig.  45.   Coordinate  system  for  oscillating  cylinder. 
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A  comparison  of  Eqs.  (46)  and  (15)  shows  that  z  is  replaced  by  the  rela- 
tive position  z  -z   in  the  latter  equation  and  two  new  terms  are  added 
v  n  o  M 

which  account  for  the  forces  arising  from  the  velocity  and  acceleration 
of  the  cy I i  nder. 

The  calculations  for  the  boundary  layer  and  the  separation  points  re- 
quire the  evaluation  of  the  tangential  velocities  relative  to  the  cylinder. 
At  an  arbitrary  time  t  and  angle  9  (measured  as  shown  in  Fig.  46),  the 
tangential  velocity  is  given  by 


o  j-r  i  lQ\        :  i  -i(8-Tr/2), 

u  =  Re  il_q(z  +  e   )  -  z  J  e        } 


(47) 
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Fig.  46.   Relative  tangential  velocity  on  boundary  of 
movi  ng  cy I i  nder. 


As  before,  the  velocity  at  any  point  in  the  flow  field  is  given  by 
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C.   NUMERICAL  PROCEDURE 

It  is  apparent  from  the  formulation  of  the  problem  that  the  solution 
of  Eq.  (44)  requires  two  parameters  (£  and  a  )  characterizing  the  mechani- 
cal system  and  the  instantaneous  fluid  force  imparted  to  the  cylinder. 
The  parameter  a  ,  known  as  the  mass  parameter,  may  be  readily  calculated 
for  a  given  cylinder  of  mass  m,  radius  c,  and  fluid  density  p,  Usee  Eq. 
C43)~J.  The  damping  ratio,  £,  is  a  measure  of  the  relative  energy  loss 
within  the  material  due  to  molecular  motion  and  does  not  include  the  so 
called  fluid  damping.   All  fluid  forces  acting  on  the  cylinder  are  included 
in  the  proper  specification  of  C,  which  represents  the  integrated  effect 
of  the  surface-pressure  loading  by  the  fluid  on  the  vibrating  cylinder. 

The  discrete  vortex  model  developed  for  the  stationary  cylinder  is 
used  together  with  Eqs.  (46)  and  (48)  to  determine  the  fluid  forces  acting 
on  the  cylinder.   Equation  (44)  is  integrated  through  the  use  of  a  fourth 
order  Runge-Kutta  method.  Otherwise,  no  additional  modifications  are 
introduced  into  the  procedures  described  in  the  previous  chapter.   In  fact, 
the  only  specification  made  was  the  time  at  which  the  cylinder  was  allowed 
to  respond  to  the  fluid  forces.   Such  a  decision  was  necessary  to  avoid 
the  coupling  of  two  transient  states  resulting  from  the  impulsive  start 
of  the  flow  and  the  transient  response  of  the  cylinder.   This  would  have 
made  the  calculations  unnecessarily  complex.   Instead,  it  was  decided 
first  to  establish  a  flow  with  a  sufficient  number  of  vortices  in  the 
wake  (e.g.,  by  regarding  the  cylinder  as  stationary  for  t  <  10)  and  there- 
by studying  the  response  of  the  cylinder  through  a  more  or  less  steady 
uniform  flow.   The  approximations  signified  by  "more  or  less"  are  due  to 
the  fact  in  all  numerical  ca I cu I  at  ions  the  vortex  street  is  finite.   In  the 
present  calculations,  the  cylinder  is  allowed  to  respond  after  t  >  10, 
which  was  sufficiently  large  for  all  intents  and  purposes. 
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Mathematically  speaking,  there  are  no  particular  restrictions  on  the 
selection  of  the  numerical  values  of  5  and  a  other  than  the  fact  that 
they  should  be  physically  realizable.  From  the  point  of  view  of  actual 
calculation,  however,  one  must  place  some  restrictions  on  £  and  a  to 
limit  the  computation  time.  The  smaller  the  value  of  5,  the  larger  the 
amplitude  of  oscillation.  The  larger  the  mass  of  the  cylinder  or  the 
smaller  the  parameter  a   (keeping  everything  else  constant),  the  larger 
is  the  transient  time  during  which  the  cylinder  oscillations  asymptoti- 
cally approach  their  final  values.   Consequently,  a  very  small  value  of 
5  coupled  with  a  relatively  small  value  of  a  would  require  many  cycles 
of  oscillation  before  the  cylinder  reaches  its  final  amplitude.   This 
makes  the  computation  prohibitively  expensive  and  time  consuming.   The 
foregoing  arguments  constituted  one  of  the  primary  conditions  in  the  selec- 
tion of  X,   and  a  .   The  second  and  probably  more  important  reason  from  a 
scientific  point  of  view  was  the  focusing  of  attention  on  the  understand- 
ing of  the  physics  of  the  phenomenon  rather  than  on  the  solution  of 
specific  problems.   With  the  foregoing  considerations  in  mind,  the  £  and 

a  values  shown  in  Table  II  were  used  in  the  calculations.   The  role  of 
o 

the  new  parameter  Sn  =   £/a  appearing  in  the  table  is  discussed  later. 

b        O 


TABLE  I  I 


SET 

? 

a 
0 

SG 

1 

0.02 

0.0125 

1.6 

1  1 

0.04 

0.0250 

1  .6 

1  1  1 

0.08 

0.0500 

1.6 

IV 

0.04 

0.0125 

3.2 
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D.   DISCUSSION  OF  RESULTS 
I .   I ntroduction 

The  major  part  of  the  discussion  of  the  results  deals  with  the 
phenomenologica I  predictions  of  the  model  and  their  comparison  with  those 
observed  or  measured  experimentally.   For  this  purpose  only  one  set  of 
damping  ratio  and  mass  parameter  is  considered.  The  variations  of  the 
lift  coefficient,  relative  amplitude,  phase  angle  between  the  lift  force 
and  the  cy I inder  response,  and  final ly  the  ratio  of  the  cy I inder-osci I la- 
tion  frequency  to  natural  frequency  is  presented  as  a  function  of  U 
defined  by 


U  = 


r   2ir  S 


(49) 


In  dimensional  terms,  U  =  U/oo  D  .  The  reasons  for  the  choice  of  this 

'   r      n 

particular  parameter  are  twofold.   Firstly,  it  is  commonly  referred  to 

in  the  previous  works  [15].      Secondly,  it  is  a  measure  of  the  velocity 

of  flow  past  a  cylinder  of  given  diameter  and  natural  frequency. 

Figure  (47)  shows  the  results  obtained  by  Feng  [11']   with  a  lightly 

damped  circular  cylinder.   As  noted  earlier,  when  the  vortex  shedding 

frequency,  f   ,  for  a  stationary  cylinder  approaches  the  natural  frequency 

o 
Cu  =  0.3  in  Fig.  (47)],  the  vortex  shedding  frequency  becomes  nearly 

identical  to  the  natural  frequency  of  the  cylinder.   Furthermore,  the 

frequency  of  oscillation  of  the  cylinder,  oj  ,  nearly  coincides  with  its 

natural  frequency.  As  U   increases  towards  1.0,  the  amplitude  of  both 

the  oscillation  and  the  transverse  force  increases  rapidly  while  u/u)  and 

n 

U  /u)  remain  nearly  constant  at  a  value  slightly  under  unity.   As  U   in- 
c  n  '  =   /  i  r 

creases  further,  both  amplitudes  drop  either  gradually  or  abruptly  as 
illustrated  in  Fig.  (47).   The  abrupt  drop  is  a  consequence  of  hysteresis 
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Fig.  47.   Experimental  results  for  self-excited  transverse  oscilla- 
tions of  a  cylinder  C77]. 
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noted  earlier.   For  U  >  I.I  the  frequency  of  vortex  shedding,  to  ,  jumps 

to  its  stationary-cylinder  value,  i.e.,  oj  =  to     However,  the  cylinder 

o 
continues  to  oscillate  at  the  natural  frequency  of  the  cylinder  with  very 

small  amplitudes.  One  can,  therefore,  identify  two  major  regions  in  the 

relationship  between  to  ,  w  ,  and  to  .   For  U  sma I  I er  than  about  0.8, 
K         v'  c'      n       r  ' 


O)      = 


U)       and  to     -  w  .      In  the   range  0.8  <   U     <    I  .  I ,   to     <  oj       and  oj     *  to     =  to   . 
v  en  3  ~r~'vv  ven 

o  o 

For    I.I    <   U     <    1.3,  to     =  to       and  to     -  to   . 
~r~  .  v  v  c  n 

o 

The  phase  angle  between  the  exciting  force  and  the  cylinder  re- 
sponse undergoes  dramatic  changes  in  the  region  where  0.8  <  U  <  I.I. 
The  reasons  leading  to  the  amplification  of  the  oscillations  of  the 
cylinder  and  the  lift  force,  and  the  rapid  increase  of  the  phase  angle 
from  about  20  degrees  to  100  degrees  have  not  been  understood  because  of 
the  difficulty  of  measuring  or  observing  the  detailed  kinematics  of  the 
flow  field  about  the  cylinder.   Many  experiments  have  essentially  confirmed 
the  general  features  of  the  results  obtained  by  Feng  [77].   Most  of  the 
attention  has  been  devoted  to  the  understanding  of  the  locking-on  of  the 
vortex  shedding  frequency  to  the  natural  frequency  of  the  cylinder  in  the 
range  of  U  values  0.8  <  U  <  I.I.   In  conformity  with  the  previous 
literature,  this  region  is  variously  referred  to  as  the  synchronization 
region,  locking-on  region,  capture  region,  or  simply  as  the  resonant  region. 

The  numerical  results  obtained  with  c   =  0.02  and  a  =  0.0125  are 

o 

shown  in  Fig.  (48),  again  as  -a  function  of  U  .   It  is  immediately  apparent 

that  the  results  shown  in  this  figure  follow  essentially  the  same  trends 

as  those  in  Fig.  (47).   Also  shown  in  Fig.  (48)  are  the  corresponding 

values  of  oj  and  ratios  to  /oj*  or  U  /U*  where  oj*  and  U*  are  the  values  of 
o  o  o     r  r       or 


' 


to  and  U  when  n.,  reaches  its  maximum  value  (point  A).   This  point  is 
referred  to  as  the  point  of  perfect  synchronization. 
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Fig.  43.   Numerical  results  for  self-excited  transverse 

oscillations  with  £=0.02  and  a  =0.0125. 
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2.  Mechanism  of  Synchronization 

Previous  investigators  L~75H  have  uncovered  a  number  of  relatively 
unrelated  facts  regarding  the  kinematic  and  dynamic  behavior  of  flow 
about  a  transversely-oscillating  cylinder.  The  most  notable  of  the  ob- 
servations and  measurements  may  be  summarized  as  follows: 
Ci)     Oscillations  increase  the  spanwise  correlation; 
Cii)    The  strength  of  the  vortices  increases; 
(iii)   The  mean  base  pressure  decreases  and  the  mean  drag  coefficient 

increases; 
(iv)    Transverse  force  is  amplified  relative  to  its  stationary-cylinder 

va I ue; 
(v)     The  phase  angle  between  lift  and  cylinder  motion  increases  with 


U  ; 
r 


Cvi) 


The  velocity  range  over  which  synchronization  occurs  increases 
with  oscillation  amplitude; 

(vii)   Separation  point  excursion  increases  with  oscillation; 

(viii)   Lateral  spacing  of  the  vortices  decreases  with  increasing  ampli- 
tude and  is  unaffected  by  changes  in  frequency  of  oscillation; 
Wake-formation  length  decreases  systematically  with  increasing 


(Ix) 


(x) 


amplitude  of  vibration  for  U  <  U*  and  increases  for  U  >  U*; 

r          r  r          v 

The  wake  contracts  for  U  <  U*  and  expands  for  U  >  U*. 

r          r  v                        v          r 


The  facts  summarized  above  have  been  distilled  from  a  number  of 
investigations,  not  all  of  which  have  been  conducted  at  identical  Reynolds 
numbers  or  under  similar  experimental  conditions.  Consequently,  there 
are  certain  features  of  the  flow  which  have  been  observed  to  increase  in 
one  investigation  and  decrease  in  another.   Whenever  such  conflicts  arose, 
heavier  weight  was  subjectively  attached  to  the  observation  more  frequently 
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reported.  Perhaps  the  most  interesting  feature  of  all  these  isolated 
facts  is  that  it  has  not  yet  been  possible  to  establish  the  cause  and 
effect  relationship  among  them  through  a  phenomonologica I ly  and  mathe- 
matically defensible  mechanism.   In  the  following,  the  mechanism  suggested 
by  the  numerical  experiments  is  described  in  detail. 

Consider  an  elastically  mounted  circular  cylinder  which  is  at  its 
lowest  displacement  during  a  given  cycle.  Also  assume  that  the  stationary 
cylinder  value  of  the  vortex  shedding  frequency  is  slightly  greater  than 
the  natural  frequency  of  the  cylinder.   As  soon  as  the  cylinder  begins 
to  move  upwards,  the  lower  attached-vortex  sheet  lags  behind  the  cylinder 
because  of  its  inertia.   In  other  words,  the  wake  axis  passing  through 
the  instantaneous  center  of  the  cylinder  rotates  clockwise  (see  Fig.  (49) 
for  t  =  189).   In  the  meantime,  the  velocity  relative  to  the  cylinder 
rotates  clockwise  an  angle  dictated  by  the  velocity  of  the  uniform  flow 
and  the  instantaneous  velocity  of  the  cylinder.   The  rotation  of  the 
velocity  vector  coupled  with  the  rotation  of  the  wake  axis  moves  the 
upper  separation  point  further  downstream  and  the  lower  separation  point 
further  upstream  relative  to  that  which  would  have  existed  on  a  stationary 
cylinder.   The  rate  of  circulation  at  the  upper  separation  point  increases 
rapidly  and  decreases  at  the  lower  separation  point.   This  gives  rise  to 
a  larger  circulation  per  unit  time  than  that  for  the  stationary  cylinder 
for  the  vortex  towards  which  the  cylinder  is  moving.   Furthermore,  the 
vortex  spiral  which  is  attached  to  the  upper  separation  point  does  not 
move  out  of  the  cylinder's  way  because  of  its  inertia  and  consequently 
comes  closer  to  the  approaching  cylinder.   These  two  facts,  namely  the 
increase  of  circulation  and  the  decrease  of  the  relative  distance  between 
the  cylinder  and  the  upper  vortex  spiral,  confirm  the  observation  made  by 
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Fig.  49.   Wake  configuration  for  self-excited  transverse  oscillations. 
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Fig.  49  (con't. ) . 


Wake  configuration  for  self-excited  transverse  oscilla- 
tions . 
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Davies  [JQl   that  "The  growing  vortex  on  vibrating  cylinders  seems  to  roll 

up  more  quickly  and  with  a  strength  about  35  percent  larger  than  that  for 

a  stationary  cylinder."  The  increase  of  the  vortex  strength  predicted 

by  the  current  model  will  be  discussed  later.   The  phenomenon  is,  however, 

relatively  more  involved  than  the  simplified  introductory  remarks  would 

suggest.  There  is  a  third  factor  which  enhances  the  strengths  of  the 

attached  vortex  sheets.  This  is  directly  related  to  the  difference  between 

the  period  of  oscillation  of  the  cylinder,  T  =  l/f  a  l/f  ,  and  the  time 
K  '      '  c     c     n 

required  for  the  shedding  of  a  vortex  from  a  stationary  cylinder,  i.e., 

T  -  T  .   Each  growing  vortex  about  an  oscillating  cylinder  continues  to 

o 
be  fed  circulation  not  only  at  a  higher  rate  but  also  for  a  longer  time 

period.   This  further  reinforces  the  strength  of  each  growing  vortex. 

It  is  now  necessary  to  establish  the  link  between  the  enhanced 

vortex  strengths  and  the  locking-on  of  the  vortex  shedding  frequency  to 

the  frequency  of  cylinder  oscillation.   The  results  of  the  stationary 

cylinder  have  shown  that  the  larger  the  strength  of  the  vortices  in  the 

near  wake  the  smaller  their  velocity  relative  to  the  cylinder.   That  is, 

strong  vortices  with  their  equally  strong  images  tend  to  conserve  their 

moment  of  circulation.   Consequently,  the  vortices  are  shed  at  larger 

time  intervals.   This,  in  turn,  leads  to  smaller  Strouhal  numbers.   The 

decrease  of  S,  or  the  increase  of  T  ,  brinqs  T  closer  to  T  .   This  ten- 

v      3   v  c 

dency  is  reinforced  every  succeedinq  cycle  until  T  =  T  -  T  or  until 

'  a   /  v    c    n 

the  effects  which  cause  the  changes  in  amplitude,  period,  lift  force, 
and  the  phase  angle  disappear.  At  such  time,  a  stable  periodic  state  is 
reached. 

The  foregoing  explanation  dealt  only  with  the  kinematic  aspects 
of  the  fluid  motion.   In  the  final  analysis,  it  is  the  balance  between 
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the  work  done  by  the  fluid  on  the  cylinder  and  the  work  done  by  the  dissi- 

pative  forces  (internal  friction)  which  determine  the  stability  of  the 

oscillations.   If  the  net  work  becomes  negative,  then  the  oscillations 

decay  and  vice  versa.  The  parameter  which  establishes  the  relationship 

between  the  friction  mechanism  and  the  given  state  of  oscillation  is  the 

phase  angle  between  the  exciting  force  and  displacement.   It  is  given  by 

[79] 

2c,   co/co 

.   tan  <j>  =  2__  (50) 

I -(co/co  ) 
n 

As  co/co  approaches  unity  from  co/co  <  I  ,  cj>  increases  rapidly,  the  rate  of 

increase  being  a  strong  function  of  £.   For  a  cylinder  with  very  small 

internal  damping,  <j>  increases  from  say  20  degrees  to  <j>  =  150  degrees  for 

a  very  sma  I  I  change  in  co/co  . 
'  3        n 

The  relationship  between  the  dynamics  and  the  kinematics  of  the 
motion  for  the  three  states  of  oscillation  and  hence  the  reason  for  the 
sustained  oscillations  in  the  region  of  synchronization  may  be  explained 
as  follows.   In  the  case  of  synchronization  (see  Fig.  50)  the  maximum 
velocity  of  the  cylinder  leads  the  maximum  lift  and  the  lift  reaches  zero 
at  about  the  time  the  cylinder  reaches  its  peak  amplitude.   In  other 
words,  the  lift  force  opposes  the  cylinder  motion  for  a  short  time  period 
as  the  amplitude  of  the  cylinder  begins  to  decrease  from  its  maximum 
displacement.   If  for  any  accidental  reason  the  amplitude  of  fhe  cylinder 
in  one  direction  were  to  be  increased  while  the  cylinder  was  undergoing 
synchronized  oscillations,  the  increased  strength  of  the  vortex  growing 
ahead  of  the  cylinder  motion  increases  the  period,  T  .   This,  in  turn, 
decreases  the  phase  angle  between  the  force  and  the  displacement  which 
then  increases  the  phase  angle  between  the  force  and  the  velocity. 
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Fig.  50.   Lift  coefficient  and  cylinder  velocity  versus  time. 
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Consequently,  the  energy  transferred  from  the  fluid  to  the  cylinder  de- 
creases below  that  dissipated  by  the  material  damping.  These  sequences 
of  events  decrease  the  amplitude  to  the  level  where  the  energy  extracted 
from  the  fluid  over  a  cycle  is  just  balanced  by  that  dissipated  by 
internal  friction. 

The  role  played  by  all  of  the  terms  of  the  governing  equations  of 
motion  and  the  energy  transfer  mechanism  may  be  quantified  as  follows. 
Re-writing  Eq.  (44) 


2ttcS  ttS      n  o     o 

n  +  n  +   ( )      11  =  217  5     a   C, 

oj  oj  o     o   L 

o  o 


(44   repeat) 


and  multiplying  with  n  and  integrating  over  a  cycle,  one  has 


t+T  t+T 

2tt<;S 


t+T 


/ 


/ 


t 


ITS 


nrid  ; °     I      f,2(jt  +( — 2.)2  I      r^dt  =   2tt2S   2a 

u)  /  a)  I  o     o 

o       J  o 


i  ™ 


t+T 


h 


ndt 


t 


(51  ) 


which  reduces  to 


t+T 


2ircS 


CO 


/ 


•2  2     2 

n  dt  =  2tt  S     a 
o     o 


t+T 


(52) 


In  the  case  of  a  perfect  balance  between  the  energy  input  Cright  hand  side 
of  Eq.  (52)^  and  the  energy  dissipated  (left  hand  side),  the  above  equality 
holds  true.   For  an  unstable  system  where  the  equality  does  not  hold  true, 


t+T 


t+T 


f    r   'h+   27T'S°  f    '< 


E  =  2tt  S  a  /  C.ndt /  n  dt 

o  o  /    L       co. 


(53) 
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Thus,  E  =  0.0  corresponds  to  steady  state  with  constant  amplitude  of  os- 
cillation; E  >  0.0  corresponds  to  amplification;  and  finally,  E  <  0.0 
corresponds  to  damping  of  oscillations.  The  variation  of  n,  n,,  v\   and  C. 
in  Eq.  (44)  as  a  function  of  time  is  shown  in  Figs.  5la-c.   Figure  51a 
corresponds  to  the  case  where  E  >  0.0;  Fig.  51b  to  E  =  0.0;  and  Fig.  51c 
to  E  <  0.0.  Clearly,  the  phase  angle  between  the  lift  force  and  the 
velocity  is  smaller  for  the  case  of  E  >  0.0  than  for  E  =  0.0.   In  fact, 
one  can  state  that  the  said  phase  angle  increases  from  a  plus  to  a  minus 
value.  Only  in  the  case  of  Fig.  51b  is  the  phase  angle  such  that  the 
energy  imbalance  is  for  all  intents  and  purposes  zero.   Furthermore,  the 
stability  of  the  system  to  small  perturbations  in  amplitude  is  such  that 
a  small  increase  in  amplitude  leads  to  a  decrease  in  phase  between  lift 
and  displacement  and  thus  to  an  increase  in  phase  between  lift  and 
velocity.   This,  as  noted  earlier,  leads  to  a  decrease  in  amplitude.   The 
reverse  is  true  for  a  slight  decrease  in  amplitude  from  a  steady  state 
osci I  I ation. 

3.   Parametric  Study  of  Synchronization 

Numerical  experiments  have  been  performed  for  the  c,   and  a  combi- 
nations given  in  Table  II  in  addition  to  that  already  presented.   The 
results  for  each  set  are  presented  in  Figs.  52  through  54.   All  of  the 
parameters  shown  in  these  figures  display  similar  variations  with  oj  or 
U  .   Furthermore,  they  are  quite  similar,  at  least  in  general  form,  to 
those  experimentally  obtained  by  Feng  (see  Fig.  47).  A  closer  look  at 
the  experimental  and  calculated  values  show  some  differences,  part  of 
which  may  be  attributable  to  the  limitations  of  the  numerical  model  and 
part  to  the  limitations  of  the  experimentation.   In  other  words,  all  of 
the  differences  to  be  pointed  out  are  not  entirely  consequences  of  the 
numer ica I  model . 
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• 

Fig.  51a.   Variation  of  n,  n,  H  and  C,  versus  time  for  E>0.0. 
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Fig.  51b.   Variation  of  n. ,  1»n  and  C  versus  time  for  E=0.0. 
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Fig.  51c.   Variation  of  r\,   n,  fi,   and  C  versus  time  for  E<0.0. 


145 


J 

0.6 
Fig.  52. 


0.7 


0.8 


0.9 


Numerical  results  for  self-excited  transverse 

oscillations  with  £=0.04  and  a  =0.025. 

o 


!46 


I 

0.6 

Fig.  53. 


0.7 


0.3 


0.9 


I  .0 


Numerical  results  for  self-excited  transverse 

oscillations  with  5=0.08  and  a  =0.05. 

o 


147 


200  p 

160  - 

120 

80 

40 
0 


0.9 

0.8  - 
0.7 

0.6 
0.5 

0.4 
0.3 
0.2 
0.  I 


0.0 


0.8 


0.9 


I  .0 


I  .  I 


.2 


0) 

o 

J l 

I  .3 


0.6 


0.7 


0.9 


.0 


Fig.  54.  Numerical  results  for  self-excited  transverse 

oscillations  with  r=0.04  and  a  =0.0125. 

*  o 


148 


A  comparison  of  Figs.  47  and  52  shows  that  the  numerical  model 
predicts  a  peak  lift  amplitude  at  a  U  value  smaller  than  U*.   In  Feng's 
experiments,  the  lift  and  oscillation  amplitudes  reached  their  peak  at 
nearly  identical  U  values.  The  reasons  for  this  are  not  entirely  clear. 
It  may  be  conjectured  that  it  is  a  consequence  of  the  extremely  low  damp- 
ing encountered  in  Feng's  experiments  (5  -  0.001).   In  fact,  mechanical 
systems  with  very  small  damping  are  much  more  finely  tuned  and  more  pre- 
cariously balanced  than  those  with  larger  internal  damping.   It  is  there- 
fore natural  to  expect  that  the  more  stable  systems  experimented  with  in 
the  present  numerical  calculations  should  exhibit  larger  differences  in 
U  values  at  which  they  reach  their  extremals. 

The  ratio  of  the  cylinder  oscillation  frequency  to  the  natural 

CO      CO 
C      V 

frequency,  —  -  — ,  remains  practically  constant  at  a  value  slightly  under 

n    n 
unity  in  the  experiments  reported  by  Feng  and  others.   In  the  numerical 

results,  co  /co  gradually  increases  from  about  0.95  to  0.99  and  then 
'  c  n  a       ' 

rapidly  increases  toward  co  /co  .  The  differences  between  the  measurements 
r   '  v   n 

o 

and  calculations  may  be  due  to  the  difficulty  of  modelling  some  secondary 

features  of   flow  which   adjust  themselves   to   render  co  /co     practically 

J  en1' 

constant. 

The  amplitudes  of  the  lift  coefficient  and  cylinder  response  cannot 
be  compared  with  the  experimental  data  because  of  the  large  differences 
in  damping  and  mass  parameters.   However,  it  is  interesting  to  note  that 
the  rati  o  of  the  maximum  I i  ft  coef f  icient  for  the  osci  I  I ati  ng  cy I i  nders 
to  the  lift  coefficient  for  the  stationary  cylinder  (standard  run)  varies 
from  1.3  to  1.6  (see  Figs.  48  and  52-54).   This  lift  amplification  is  in 
conformity  with  most  of  the  data  reported  in  the  literature  C"75].   The 
cy  I  inder  response  shows  a  dramatic  drop  for  a  sma  I  I  change  i  n  co  or  U  , 
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reminiscent  of  the  hysteritic  drop  in  Feng's  data.   It  should,  however,  be 
emphasized  that  the  drop  inn  in  the  numerical  calculations  has  nothing  to 
do  with  the  hysteresis  resulting  from  the  change  in  damping  of  the  cylin- 
der.  The  sudden  drop  is  a  consequence  of  the  changes  in  behavior  of 
flow  external  to  the  cylinder.  One  further  distinction  must  be  made  between 
the  experimentally  obtained  sudden  drop  and  that  obtained  numerically.   In 
the  experiments,  the  velocity  of  the  flow  (or  U  )  is  increased  gradually 
through  the  U  value  at  which  the  sudden  change  occurs.   In  the  numerical 
calculations,  U  or  oj  is  kept  constant  at  each  point.  Consequently,  what 
is  observed  is  not  a  sudden  change  in  time  but  a  change  in  amplitude  as 

a  function  of  to  . 

o 

The  cylinder  response  for  oj  >  go*  does  not  remain  at  a  constant 
'         r         o    o 

amplitude.   It  is  because  of  this  fact  that  part  of  the  response  curves 

corresponding  to  that  region  are  shown  by  broken  lines  in  Figs.  48  and 

52-54.   A  sample  plot  of  n  for  to  >  oj*  i  s  shown  in  Fig.  55. 
r   K  o    o  a 

In  discussing  the  material  damping  and  mass  parameters,  another 
parameter  denoted  by  S„  =  £/a   has  been  introduced.   Among  others,  Vickery 
and  Watkins  C8CG  have  attempted  to  show  that  r\*    (r\     at  perfect  synchroni- 
zation) is  uniquely  determined  by  the  value  of  a  response  or  stability 
parameter,  S„.   Sarpkaya  l8  I J  has  shown  that  n.*,  may  be  fairly  accurately 

b  M 

represented   by 

cL   (.*)2 

njj--£— —  (54) 

/0.06   +  S„2 

b 


where  C   is  the  lift  coefficient  for  the  stationary  cylinder.   The  values 

o 

of  5  and  a  were  judiciously  selected  to  arrive  at  identical  values  of  S„ 
O  b 

for  three  sets  of  calculations  and  twice  the  S„  value  for  one  set  of 
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calculations.  According  to  the  hypothesis  advanced  by  Vickery  and  Watkins, 
n*  should  be  uniquely  determinable  by  S„.  The  results  obtained  through 
use  of  the  model  are  tabulated  below.  The  results  in  Table  III  show  that 
cylinder  response  at  perfect  synchronization  is  not  very  sensitive  to  the 
individual  variations  of  the  damping  and  mass  ratios. 

TABLE  I  I  I 


SET 

5 

a 
o 

SG 

*         

From  Eq.  (54) 

Numerical  Results 

Eq.  (56) 

1 

0.02 

0.0125 

1.6 

0.52 

0.54 

0.56 

1  1 

0.04 

0.0250 

1.6 

0.58 

0.58 

0.58 

1  1  1 

0.08 

0.050 

1.6 

0.52 

0.61 

0.58 

IV 

0.04 

0.0125 

3.2 

0.24 

0.25 

0.26 

An  exact  relationship  between  S_  and  other  parameters  may  be 
established  through  the  use  of  the  energy  balance  at  perfect  synchroniza- 
tion.  Re-writing  Eq.  (52) 


t+T 


t+T 


27TCS 


Ui 


!■ 


2  2   2 

n  dt  =  2ti ■  S  a 

o  o 


/   cLndt 


(52  repeat) 


and   expressing   £  and   C.    as 


n  =  njjj  sin  oo*t   ,   CL  =  CL*  sin    (u*t  +  <j>*) 


where  $*    is  the  phase  angle  between  the  force  and  the  displacement  at 
perfect  synchronization  one  has 
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t+T 

lidt 


/  ^ 


Trcj*S   C.  *si  nd>* 

r  £  *c        ">  _OOL  ,«, 


/** 


M 


or 


7tw*S  C.  *s  1  neb* 

n*  = 5?c (56) 

b 


Equation  (56)  shows  that  n^  is  not  a  unique  function  of  S_.   However, 
since  cb*,  C,  *,  and  co*/cj*  do  not  vary  significantly,  n,y  varies  approxi- 
mately as  S„   .   A  comparison  of  Eqns.  (54)  and  (56)  shows  that  Eq .  (56) 

b 

requires  the  use  of  parameters  which  must  be  determined  either  numerically 

or  experimentally  whereas  Eq.  (54)  requires  only  a  knowledge  of  C.   and 

o 

oj  or  U  . 
o     r 

The  results  obtained  through  the  use  of  Eq.  (56)  for  the  four 
sets  of  input  parameters  are  also  shown  in  Table  III.   Clearly  the  direct 
and  indirect  predictions  of  n*  are  in  excellent  agreement,  as  would  be 
expected. 

The  n^-Sr  relationship  given  by  Eq.  (56)  is  expected  to  hold  true 

for  values  of  co  other  than  uj*  provided  to   is  confined  to  the  region  of 
o  00 

synchronization.   Any  small  deviations  between  the  predictions  of 


ttci)  S  C,  si  neb 

o  o  L   r  ,  c,  > 

n.„  =  c (56a) 

b 


from  those  directly  calculated  with  the  model  are  simply  due  to  small 
modulations  in  the  amplitude  of  oscillations.   Numerical  calculations  which 
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can  easily  be  performed  through  the  use  of  Figs.  48  and  52-54  have  shown 
that  the  predictions  of  Eq.(56a)  are  within  a  few  percent  of  those  obtained 
di  rectly. 

The  transverse  oscillations  of  a  cylinder  in  uniform  flow  increase 
the  affective  projected  area  of  the  cylinder  and  hence  the  bluffness  of 
the  body.  Consequently,  hydroelastic  oscillations  of  a  cylinder,  such  as 
the  strumming  of  a  cable,  not  only  cause  alternating  stress,  fatigue, 
and  noise  but  also  increased  drag. 

The  numerical  model  calculated  two  important  characteristics  of 

the  fluctuating  drag  force:  the  mean  drag  and  the  amplitude  of  oscillation. 

Figure  56  shows  the  mean  drag  coefficient,  CR,  and  the  amplitude  of  the 

fluctuating   component,    ACL,    as   a   function   of  to     or   U    .      Also   shown    in 
3  r  D  or 

Fig.  56  is  a  data  point  obtained  experimentally  £8  I  J, for  the  corresponding 

to  and  n.,,  va  I  ues. 
o      M 

Figure  57  i s  a  representative  plot  of  the  instantaneous  value  of 
C  as  a  function  of  time.   It  shows,  among  other  things,  that  the  drag 
force  oscillates  at  a  frequency  twice  that  of  the  lift  force  (for  lift 
T  =  10  and  for  drag  T  ~   5). 

The  separation  points  on  an  oscillating  cylinder  oscillate  with 
a  frequency  equal  to  that  of  the  cylinder.   Representative  results  pre- 
sented in  Fig.  58  show  that  the  mean  position  of  the  separation  angle 
remains  practical ly  constant  at  a  value  of  8  =  77  degrees.   The  amp  I  itude 
of  the  fluctuations  follow  a  pattern  very  similar  to  that  of  n,M-   The 
stagnation  point  undergoes  s imi lar  osci  I  I  at  ions.   The  results  have  shown 
that  the  instantaneous  position  of  the  stagnation  point  nearly  coincides 
with  that  which  would  have  resulted  with  a  uniform  flow  at  the  instantaneous 
total  velocity.   In  other  words,  the  small  oscillations  of  the  stagnation 
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Fig.  56.   Drag  coefficient  versus  cj  for  £=0.02,  a  =0.0  f25 . 
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point  relative  to  the  uniform  flow  are  obscured  by  the  larger  oscillations 
resulting  from  the  oscillations  of  the  cylinder  or  of  the  ambient  flow 
relative  to  the  cylinder. 
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VI  I  .  CONCLUSIONS 

The  results  presented  herein  warrant  the  following  conclusions: 

1.  A  discrete-vortex  model  based  on  the  redi scretization  of  shear  layers 
and  the  use  of  appropriate  vorticity  cancellation  mechanisms  can  produce 
results  which  are  in  conformity  with  those  obtained  experimentally; 

2.  The  numerical  model  allows  experimentation  with  the  variation  of  any 
one  of  the  independent  parameters  and  permits  one  to  understand  the 
governing  physical  mechanisms; 

3.  It  has  been  shown  that  two  major  flow  parameters,  namely  the  drag  co- 
efficient and  the  Strouhal  number,  are  least  sensitive  to  changes  in  the 
characteristics  of  the  wake.   The  Strouhal  number  is  intrinsic  to  the 
shape  of  the  bluff  body  and  its  magnitude  cannot  be  varied  more  than  a 
few  percent  by  varying  the  strength  of  the  vortices  by  physical  means. 
Experiments  over  the  past  50  years  are  in  conformity  with  these  conclu- 
sions; 

4.  The  lift  coefficient  is  the  one  most  sensitive  to  the  changes  of 
flow  in  the  near  wake.   This  finding  of  the  numerical  model  and  the 
difficulty  of  obtaining  consistent  laboratory  data  for  lift  are  in  con- 
formity with  each  other.   The  numerical  model  explains  the  physical 
mechanism  as  to  why  lift  should  be  a  strong  function  of  the  evolution 
of  the  near  wake; 

5.  The  discrete  vortex  model,  as  all  other  numerical  methods  attempting 
to  model  a  turbulent  region  of  a  flow  and  its  interaction  with  a  laminar 
region,  is  in  need  of  a  hypothesis  to  properly  simulate  the  circulation 
dissipation  and  the  vorticity  diffusion  caused  by  turbulence. 
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Notwithstanding  this  fact,  a  very  simple  decay  mechanism  is  capable  of 
simulating  the  said  consequences  of  turbulence; 

6.  The  predictive  powers  of  the  numerical  model  have  been  put  to  test 
through  its  use  in  the  study  of  self-excited  oscillations  of  elastically 
mounted  circular  cylinders.  The  results  were  found  to  be  in  conformity 
with  practically  all  observations  and  measurements  made  in  recent  years 
and  led  to  the  understanding  of  the  mechanisms  governing  the  synchroni- 
zation phenomenon. 
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